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COULOMB CORRECTIONS TO THE BETHE-HEITLER CROSS 
SECTIONS FOR ELECTRON-NUCLEUS BREMSSTRAHLUNG 

SUMMARY 


The electron -nucleus bremsstrahlung cross section is computed by 
using the closed analytic expression for the Born approximation and a more 
elaborate calculation based upon a partial wave expansion. Theoretical details 
and a comparison of the results obtained by the two methods are presented. 

A correction factor based upon the two methods is obtained. This 
factor is then used to correct the Born approximation for Coulomb effects 
that have been neglected due to the use of plane waves for the electron in the 
matrix element. 

The corrected Born approximation gives relatively accurate results 
at the lower portion of the energy region (0. 1 to 2. 0 MeV) under investiga- 
tion. However, because of a lack of data at the upper portion of the energy 
region, the correction produces an overestimation of the cross section for 
higher energies. 

The correction factor may be extended to higher energies when data 
from the partial wave approach become available. The limitations at the 
present time are computer storage and long computation time resulting from 
the involved nature of the partial wave method. 


SECTION I. INTRODUCTION 


Radiation protection for man and sensitive instruments in the environ- 
ment of space involves two primary sources of radiation: high energy 
nucleons and the radiations trapped in the earth's magnetic field. The latter 
source, consisting of high energy charged particles, is often of the greatest 
importance because of its abundance, spatial distribution, and energy 
spectra. 

To determine the radiation that penetrates the spacecraft, one must 
have a knowledge of the radiation sources and understand the interaction 



of the radiation with matter. Also, among other processes, an evaluation of 
the penetration and transport of the electron must be made, which requires 
accurate knowledge of the basic physical processes involved. This paper is 
devoted to a study of the basic process involved in the loss of electron kinetic 
energy because of radiation emission in a medium. 

It is well known from classical electromagnetic theory that whenever 
a charge is accelerated as it passes through matter, radiation emission will 
occur. This radiation is known as bremsstrahlung (braking radiation) . The 
acceleration, in the case of electron-nucleus bremsstrahlung, is the result 
of deflecting the path of the electron as it is scattered by the repulsive 
Coulomb field of the atomic nucleus. 

The calculation of the electron bremsstrahlung cross section is a 
problem of long standing in theoretical physics. The exact solution of the 
nonrelativistic bremsstrahlung problem has been obtained by Sommerfeld 
[1]. Bethe and Heitler [2] first formulated the relativistic theory by using 
Dirac's electron theory and the Born approximation. (Throughout the 
literature, the synonymous terms Bethe-Heitler theory and Born approxima- 
tion for bremsstrahlung production are used. The same convention will be 
followed in this paper.) Schiff [3] developed a theory that included screening 
effects due to the atomic electrons, but this theory is only valid for high 
energies and small-to-moderate angles of deflection. Bethe and Maximon [4] 
have presented a theory, exclusive of the Born approximation, that is valid 
only in the extreme relativistic region above 20 MeV. Subsequent develop- 
ments have included various corrections to the Born approximation, valid in 
specialized energy ranges, that account for atomic screening and Coulomb 
effects not included in the Bethe-Heitler theory [5] . 

For incident electron energies below approximately 10 MeV, the 
Bethe-Heitler theory gives incorrect results for the cross section over the 
entire photon energy spectrum [5] . This inaccuracy can be attributed almost 
entirely to the use of plane waves in the Born approximation, since for low- 
to-moderate energies, the plane wave is distorted by the Coulomb field of the 
nucleus. Fpr example, in the region near 0. 1 to about 2. 0 MeV, the 
accuracy, when theory is compared with experiment, is only within ± 20 
percent [5] . For incident electron energy from 0. 5 to 1.0 MeV [6] and at 
1. 7 MeV [7] , the experimental data are higher than the theory predicts, 
while at 2.72 MeV [8] the theoretical results are higher than the experi- 
mental results. This would suggest that a "transition region" exists between 
the latter two energies. 
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Elwert [9] has developed an analytical Coulomb correction factor, 
valid up to approximately 0. 1 MeV, on the basis of a comparison between the 
nonrelativistic Sommerfeld theory and the nonrelativistic Born approximation. 
However, in the energy range from 0. 1 to 2. 0 MeV, Coulomb corrections to 
the Bethe-Heitler theory are not available in analytical form [5, 10] and are 
presently estimated by using a corrective multiplicative factor. This factor 
is the ratio between the experimental total radiation cross section and the 
calculated one using the Born approximation. Interpolation and extrapolation 
techniques are then used when experimental data are not available for specific 
energy regions and/or materials. 

Exact results using the relativistic theory for bremsstrahlung cross 
sections are nonexistent. However, Brysk, Zerby, and Penny [11] have 
attempted a formulation of the problem using exact Coulomb wave functions 
instead of the plane waves used for calculating the matrix elements in the 
Born approximation. This formulation essentially involves a partial wave 
expansion for the incident and scattered electron and for the photon. Since 
the incident electron energy dictates the number of matrix elements required 
for a solution of the problem, computer precision and storage limitations 
become limiting factors for higher energies. Therefore, this approach is 
practical only for low -to -mode rate energies. 

The Born approximation for bremsstrahlung cross sections is 
available in closed analytic form and is relatively easy to calculate on 
modern, high speed digital computers, while the partial wave approach is 
not. Therefore, it would be advantageous to develop a Coulomb correction 
factor to the Born approximation, valid in an energy range that is higher 
than the energy at which the Elwert factor is valid. This factor could be 
obtained based upon a comparison between the Born approximation and 
experimental data. However, it is desired to obtain a correction factor based 
upon two theoretical approaches. Therefore, the correction factor here is 
based upon a comparison between the relativistic Born approximation and 
the partial wave solution, that will be essentially an extension of the Elwert 
factor to higher energies. 

The primary energy range of interest is the "transition region. " 
Therefore, our efforts have been applied to the area of 0. 1 to 2. 0 MeV with 
possible extensions to higher and/or lower energies. 
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SECTION II. THEORY 


A. The General Bremsstrahlung Problem 

The bremsstrahlung problem may be formulated as follows: Assume 
that an electron of known kinetic energy collides with a thin material of 
atomic number Z. Find the probability that a photon with energy between 
given limits will be emitted into a solid angle measured with respect to the 
direction of the incident electron. A more general form of the problem would 
be to determine the polarization of the photon and scattered electron. How- 
ever, for purposes herein, it will be assumed that the incident electrons are 
unpolarized and the polarizations of the emitted or scattered particles will 
not be of interest. 

The events occurring in the process are illustrated in Figure 1. An 
electron, incident along the vector p^, is deflected by the nuclear Coulomb 

field of the scattering atom with the resulting emission of a photon along the 
direction k and the scattering of the electron along direction p . 

The problem of calculating the bremsstrahlung cross section can be 
formulated by two principal methods. First, the Born approximation can be 
utilized with the results limited to the applicability of the Born approximation; 
i.e., 27TZ/137 « /3 q , /?, where (/3) is the ratio of the incident ( scattered) 

electron velocity to the velocity of light. This restriction is relatively severe 
for low energies and/or high atomic numbers. On the other hand, one could 
attempt to calculate the cross section using exact Coulomb wave functions 
instead of the plane waves for the electron initial and final states. This is a 
more difficult approach. In fact, the scattering problem using the Dirac 
wave equation with a Coulomb field has not been solved in closed form. The 
second approach could be approximated by using partial wave expansions with 
various numerical techniques. 

To illustrate the difficulty encountered in obtaining the exact solution, 
one first assumes that the range of the force is on the order of the dimensions 
of the atom. The number of important values of angular momentum that must 
be included is then approximately b/X, where b is the Bohr radius and X 
is the de Broglie wavelength of the incident electron divided by 27r. It can be 
shown [4] for the extreme relativistic case that b/X is approximately equal 
to 137 E^, where is the incident energy in m^c 2 units. For example, 
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k (EMITTED PHOTON) 



(INCIDENT ELECTRON) 


Figure 1. Geometry of the bremsstrahlung problem. 


consider an electron with a kinetic energy of 20 MeV; this gives 


b/< » (20/0.511) 137 , 


or 5480 important values of angular momentum to be considered and 
approximately 3 x 10 7 matrix elements required for the exact solution. 


The obvious advantages and disadvantages of the two methods, in the 
energy range of interest, are that the former is simpler but less accurate, 
while the latter is more accurate but more difficult to apply. 

To present the theoretical results necessary to obtain the desired 
correction factor to the Born approximation, one first derives the expression 
for the Born approximation using Feynman's method [12] to obtain the matrix 
elements. The Coulomb wave function partial-wave expansion is then 
formulated beginning with the separation of the Dirac equation into polar 
coordinates, solving the radial equations, and then making the necessary 
expansions to obtain a calculable form for the cross section. 


UNITS, SYMBOLS, AND NOTATION 

Throughout the following development, unless otherwise specified, 

nuclear dimensionless units will be used where fi = m = c = 1. The units of 

o 

energy will then be m^c 2 ; momentum, m^c; angular momentum, 'fi; 
length, -ft/m^c; and the cross section will be in units of (ft/m^c) 2 . 

As an illustration of the units and a calculation of some required 
quantities , one notes that the total energy for a system is given by 
E 2 = m Q 2 c 4 + p 2 c 2 , where is the particle rest mass and p is the 

momentum. Denoting, for the present, by primes the energy in m c 2 
units, one writes 

e' 2 = E 2 / (m c 2 ) 2 = i + p 2 c 2 /m 2 c 4 , 

0 o 

or e' 2 = i +p 2 /(m Q c) 2 . However, the units of p are m^c; then 

e' 2 = 1 + p' 2 . Hence, the momentum can be calculated from the simple 
relationship p 2 = E' 2 - 1 . 

Relativistically, the energy is given by 
E = m c 2 = m c 2 ( 1 - j3 2 ) ^ , 

O 0 *0 
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where B =v /c ; v = the velocity of the particle. Likewise, 

p = m v ( 1 - /3 2 ) 1//2 = m c /3 ( 1 - /3 2 ) 1//2 or p ' = j8 (1-/3 2 ) 

*o o o o o o o o o o 

Thus, Pq/Eq = 0 and /3 = p/E, where P Q (p) denotes the initial 

(final) momentum, etc. (The primes have been dropped, since these units 
will be used henceforth unless otherwise specified.) 

The following symbols will be used unless otherwise indicated: 

p^Cp) = initial (final) electron momentum 

E q (E) = total initial (final) electron energy 

k ( k) = photon energy ( momentum) 

/3 (/?) = ratio of initial (final) velocity of electron to the velocity of 

° light, c 

m = rest mass of the electron, 0. 5110062 MeV 
o 

e 2 /fic = a = the fine structure constant 1/137. 0367 

O/O — 13 

r = e /m c = classical electron radius, 2.82x 10 cm 
o o 

T (T) = initial (final) electron kinetic energy 

6 Q ,9,<p are the angles defined in Figure 1 

dcr = the differential cross section for bremsstrahlung 

X = fi/m c = the Compton wavelength of the electron, 

0 ° 3. 86 x io ” 11 cm' 

The following notations will be used throughout unless otherwise 
indicated: 

r — a vector (three-dimensional) 
r — a unit vector 
a — a matrix 
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a — a matrix vector 

p — a four vector 

p — a four vector component 
M 

p^y^ = P-T^ + P 474 ) i = i, 2, 3, where y are the well-known y 
matrices; for example, see Reference 12, p. 39. 

Q* = the complex conjugate of Q 

Q + = the Hermetian conjugate of Q 

Q = the transpose conjugate of Q 

B. Cross Section for Bremsstrahiung Production 

The cross section for bremsstrahiung production can be calculated by 
using the methods of time -dependent perturbation theory. The cross section 
for one photon emission will be calculated in a cube of side L from the 
transition probability per unit time for one atom, given a current density of 
one electron per unit area per unit time in the initial state. 


do- =— dE dfi dfi, , (1) 

J y p k 

where 

da is the differential cross section for the scattering of one incident 
electron into a solid angle d £2^ and the emission of a photon, 

between the energies of E and E + d E , into the solid angle 

y y y 

dfi k 5 

J is the incident electron current density (J = |v q |/L 3 , where 
v^ is the electron incident velocity and L 3 is a volume large 
compared to other dimensions) ; 

d £2, = the element of solid angle in the direction of k . 

( d £2 , = sin 9 dd d<£ ) ; 
k o o o 
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dS2 = the element of solid angle in the direction of p . 

P ( d $2^ = sin Odd d cp) ; 

and 

w = transition probability per unit time, 
w is given by the well known expression (e.g. , see Reference 13) 



( 2 ) 


where M^. is the matrix element for the transition of the system from the 

state before photon emission to a final state occurring after photon emission 
and is the density of final states derived in, the next section. 


Equation ( 1) can be rewritten if one considers the incident velocity 

of the electron. Since B = p /E , then v = p /E , and equation ( 1) 

o o o o o o 

becomes 


E 

da = w — - L 3 d E d £2 d fl . 

P Q y p k 


After substitution of equation (2), one has 
E 

da = 2 tt — L 3 | ML. I 2 p.dE d d£2. . (3) 

p fi f y p k 

o ' 


DENSITY OF FINAL STATES 

The density of final states is defined to be the number of states 

available to a particle per unit energy interval. To calculate this, first 

consider a free particle in a cube of side L. The conditions for quantization, 

found by solving the Schrodinger equation for a particle in this cube, are 

p = n -h/L, p = n 1i/L, and p = n 'fi/L ; where p is the x component 
*x x y y z z x 

of momentum, etc. and n , n , and n are positive integers. This is 

x y z 

written as 
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p 2 = p 2 + p 2 + p 2 = ( 7 r 4 /L) 2 (n 2 +n 2 +n 2 ) 
x y z x y z 

or n 2 + n 2 + n 2 = ( pL/ 7 r 6) 2 , the equation of a sphere of radius pL/irli . 
x y z 

Now consider the octant of this sphere where all the integer 
coordinates, n., are positive. The number of the quantum states with 

momentum < p is equal to twice the volume of this octant; i.e. , 


N e =2 [(T)(f’ r ) <P L/rt > 3 ] ' < 4 > 

The factor of two is included to account for the two possible spin orientations 
of the electron. For the total number of states, the photon must also be 
included; thus, 


N = 


» 3 * s • 


( 5 ) 


Now the density of final states is 
>2 T 


9 ^N(E,E ) * t a ai 

y 8 N dp dk 


P f 9E9E 


7 


8p9k dE dE 


7 


thus, 


a 2 N 
aE a e 


- ■ V (-ikY (* £) ( 3k! I ) 

7 7 


, . dp E , dk ... 
and since = — and —— = 1 , then 

dE p dE 

7 

p { = 4 < 4 ’ r > J hszf t p2k2 • 


( 6 ) 


Since the interest here is scattering into specific solid angles df2^ and d 12^, 
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N 


instead of the integrated solid angles implied by equation (4) , one multiplies 
by (47t) . Therefore, 

^ E o E f o < 4 > lM fl l ! k 2 dk d S2p d£J k . (7-a) 

Since only unpolarized incident electrons are of interest here, and the 
polarization of the scattered electron or photon is not of concern, one writes 

T 2 E o E f “ k • (7-b) 

flfi O 

7 e 

where the factor \ comes from averaging over the polarization states of the 
incident electron and ^ implies a summation over the electron spins in 

y e 

the final state and a sum over the photon polarizations. 

C. The Bethe-Heitler Formulation of the Bremsstrahlung 

Problem 

To present an outline of the theory required to obtain the Bethe-Heitler 
formula for the bremsstrahlung differential cross section, one follows 
Feynman's [12] approach. It is not the purpose of this paper, however, to 
consider in detail the theoretical aspects underlying Feynman's procedure. 

The primary purpose is to investigate the recent calculations utilizing the 
other major approach to solving the bremsstrahlung problem; i.e. , the 
partial wave method [11] . It is hoped that this investigation will yield a 
correction factor which can be applied to the Bethe-Heitler formula and thus 
make it possible to obtain accurate results more easily than was previously 
possible. Therefore, this section is included only for completeness and the 
details not included here may be obtained from the literature [2, 12, 14] . The 
derivation of the Bethe-Heitler formula presented in this paper was, in part, 
taken from the book by Akhiezer and Berestetskii [15] . 

To illustrate Feynman's procedure for obtaining the matrix elements 
directly from a Feynman diagram, one considers the following: The matrix 
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elements in the momentum representation, for a free particle under some 
perturbation potential -ie/f(2) , are found to be [12] 

1 d* 

M = -i / u 2 e/(q 2 ) [i (pi +fii) -m] e q x ) u t . 

The various factors in the above equation are interpreted as follows 
(Fig. 2) : The particle enters the region at point 1 and moves to point 3 as a 
free particle. The particle is then scattered by a photon that has momentum 
under the action of the potential -ieX(qi). After this momentum is 
absorbed by the electron, it then moves to point 4 again as a free particle. 

A second photon scatters the electron at point 4 imparting the momentum fa. 
The particle then moves to a point 2 as a free particle with momentum 
fa + tfi + & 2 and i s described by wave function u 2 . Note that the inte- 

gral of the matrix elements is only taken over q 1 since q 2 is determined 
from the other momenta. 
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Figure 2 . Feynman diagram used to interpret the amplitude factors in the 
matrix element for a free particle-photon interaction [ 12] . 
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Feynman sets forth several rules to be followed in writing down the 
matrix element amplitude factors . 

Consider the matrix element given by 


M f . = ( u 2 Muj) . 


The rules for writing down the factors in M are [15] : 


1 . Every external electron line corresponds to a spinor 

r — *• r — 

of one of the following types u ( p) , TT ( p) , 

r r 

u and IT corresponding to the annihilation and 
creation of an electron of momentum p and polariza- 
tion r . . . 


2 . 


Every external photon line representing a photon cor- 
responds to a matrix &/ \l 2k , where k is the 
( momentum) and e is the unit polarization vector 


(X = e y 

M M 




Every external photon line repre- 


senting an external electromagnetic field corresponds 
to a matrix X(q)/(27 r) 4 . 


3. Every internal electron line of momentum p cor- 
responds to a matrix -i/ ( ipr+ m) . 


4 


5. To each vertex of the diagram there corresponds a 
6 -function containing the momenta of all the lines 
converging at this vertex, with the momenta at the 
two ends of an internal line being taken of opposite 
sign. 

6 
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8. The numerical factor which appears in M in front 

of the product of the spinors u, v, IT, v and the matrices 

y, -i/(ij^+m) is equal to ( - i) n ^<5 (2ir) ^ , where 

P 

F is the total number of internal lines, SL is the 
number of electron loops with an even number of electron 
lines. . . 

9 . Integration is carried out over the four-momenta of the 
internal lines representing virtual particles and over the 
variables q associated with the external potentials, 
and summation is carried out over the four-dimensional 
polarizations of the virtual photons. 

(The numerical factor indicated in rule 8 of Akhiezer and Berestetskii 

omitted the factor e n where e is the electron charge and n is the 
number of vertices in the Feynman diagram. ) 

The factor X(q) is proportional to the Fourier Transform of the 
potential: 

X(q) = /X(x) exp(-iq’x) dSc 

and 

A(x) = j^j-4 /X(q) exp(iq-x) d 4 q . 

For example, in the case of the Coulomb potential, one has 

A(q) = (ze/q 2 ) < 5 (q Q ) 7-4 i ( 2 -tt) 2 (2) . 

To illustrate this, first assume that one has an infinitely heavy nucleus. 

The potential is thus that of a stationary charge; i.e. , 0 = Ze/r, A = 0 , 
and X(x) = iy 4 Ze/r . Substitution of X(x) into the integral for X(q) gives 


X(q) = i ~ f exp(-iq’r) d^ f y 4 exp(iq o t) dt 
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A\q) = 1 f ~ exp(-icf- r) [y 4 <5(q Q )] ( 2 t) 

OO +i 

= ( 2w ) 2 i y 4 6 ( q ) Ze f f exp(-iqr/n) r dr d/u 
° 0-1 

oo / -iqr_ iqr \ 

= (2tt) 2 iy 4 6(q o ) Ze / (- -r^ ) r dr . 


A "convergence factor" a must now be introduced to complete the 
integration; i.e. , 


(27t) 2 iy 4 d(q Q ) Ze [ Hm 


A(q) = 


-lq 


Q!- 


,, f dr [exp( -iqr-a; ) -exp(iqr-a )] 

0 'i 


(2lr> lY 4 a(q o ) Ze | Um 
-iq | ce — ^0 


-2iq 

aT^+q 5 " 


= 2(27r) 2 iy 4 6 (q o ) Ze/q 2 


The matrix element for the cross section contains 6 -function that implies 
energy-momentum conservation. This 6 -function contains the three- 
momentum transfer to the external field, but there is no energy transfer 
since a static field has been assumed; thus, the presence of the function 
6(q^) . Therefore write V(q) , the three-dimensional Fourier transform of 

the potential, as 


V(q) =47 r Ze/q 2 , A(q)=iy 4 <5(q Q ) (27r) V(q) . 

Under the assumption that the Coulomb potential of the nucleus acts 
only once ( since this is the Born approximation) , there are two indistin- 
guishable sequences that the bremsstrahlung process may follow. The 
Feynman diagrams for these two sequences are shown in Figure 3. In the 
first case (Fig. 3a), the electron enters with momentum p^, interacts 

with the Coulomb potential of the atom, V(q) , and receives an additional 
momentum q. Later, the electron emits a photon with momentum k and is 
scattered with momentum p. The interaction responsible for the latter 
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process is denoted by e, the polarization of the photon. The second sequence 
that the process can follow (Fig. 3b) is to first emit a photon and later 
interact with the Coulomb potential. 




Figure 3. Feynman diagrams for bremsstrahlung production [12] . 

Following Feynman's approach, one defines the matrix element for 
bremsstrahlung production as M^. = (u^Mu.) or 

I M 1 2 = |(u Mu.) | 2 , 

fi f i 
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where is the matrix element given in equation ( 7-b) and ( u^) is the 

final (initial) wave function of the electron. In the Born approximation these 
are assumed to be plane waves. 

Since the final states can be obtained in two different ways, the total 
amplitude is found to be the sum of the two individual amplitudes. Thus, 

M = M( process a) + M( process b) . Using Feynman's rules, as presented 
by Akhiezer and Berestetskii [15] , one writes 


M 


and 


[luSTs] 3 02(21)1 


— TTT — — — e 2 ( 2ir ) 4 8( p +q-p-k) d 4 q 


/ 


(2tt) 4 


Therefore, 


M=M +M = 
a b 


’ * {* [ l( ** 


W) +mj 


X(q) +X(q) 


[_^ 1 




x 6(p +q-p-k) d 4 q . 
o 


The matrices in the denominators may be eliminated by noting that 
[ -iQaf+jk) + m] [ i (jgT + J^) + m] = p 2 + 2pk + m 2 = 2pk , since (iyp) 2=- p 2 , 
p 2 = -m 2 , and k 2 = 0 . Likewise, [-iQa' ~X) +m] [i(jaf "X) +m] = -2p k . 
Thus, o o o 


— f ( 

-m ' 

42k \ 

2 p k 


X(q) + X(q) 


i(y o -i^) -m 

-2 p k 
o 


& f 5(p o +q-p-k) d 4 q 


The following definitions are given to simplify our calculations [15] : 

fi = P + k = p Q + q 
f 2 = p o -k = p- q 
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Now, 


f t 2 = p 2 + 2pk + k 2 = -m 2 + 2pk + k 2 


and 


f 2 2 = -m 2 - 2p k + k 2 ; 


or since k 2 = 0 , 

f l = -'m 2 + 2pk 
and 

f 2 2 = -m 2 - 2p^k . 
One further defines 


ff + m 2 , _ f? 2 + m 2 

— and «2 = - ; 


then, 


m 2 /q = 2pk, m i! /< 2 = -2p Q k . 


, 2 .,. = 


Thus, 


M = — f l jef 

n/ 2K l 


iX - m 
m 2 /^ 


X(q) + X(q) 


iX - m 

— <s — 

m «2 J 


* ) <5(p Q +q-p-k) d 4 q 


Now X(q) = 741 V(q) 2n 6(q Q ) ; then 


-e 2 


M = (2; r) / V(q) fir 


^2k 


iX^jn 


74+74 


iX - m 

— S 

m « 2 


jef V 6(E Q -E-k) 


x <5(p o +q-p-k) d 3 q 


After integrating over d 3 q one finds 


-e 2 


M - (2tt) V(q)( d 

\[2k \ 


isfj - m 


m^c7 


74 +74 


iX ~ m" 


iiV- 

m k 2 


d J 6(E Q -E-k) 


Following Akhiezer and Berestetskii, one writes the differential cross 
section as 
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where 


da = / -j |M f .| 2 j~r 6(E Q -E-k) pEdEdS2 p d s r , 
E 


-e 


M f . = u 2 M Uj = ( 27r) V( q) u 2 Qu t 

11 *MT 


and 


Q = W 


[^r] "♦"P&r-] '| 


Then, 


Now 


then 


d<J = ^ I u 2 Q Ui 1 2 / 5(E o -E-k)|pl E dE dft p d 3 iT 


dk = k 2 dk dn 


4 . 

da = |V(q5l 2 |u 2 Q Ui 1 2 / ^(E^-E-k) |p| E dE k 2 dk d£2 p d£2^ 


k 


Integrating over E to eliminate the 6 -function gives 


e 4 — _ 1 — 

dcr = I V( q) | 2 I u 2 Q Uj 1 2 J(2?r) 5 Ip I E k 2 dk dfi dft , 

or since J = |V I = Ip l/E , the differential cross section may be 
o o o J 

written as 

4 l~^l ^ ^ 

del = ^- IV(q)l 2 IIjQu,l 2 — E— ° b k 2 dk dt! dS2 fc . ( 8-a) 

,P o' 

In computing cross sections for various processes, one first arrives 
at a cross section for particles with definite spin states. However, in our 
current problem, it is assumed that the incident electrons are unpolarized 
and the spin states of the outgoing particles are not observed. Therefore, 
one must average over the spins of the incoming electrons (since there is 
equal probability of initial spin in either direction if the electrons are 
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unpolarized) and sum over the spins of the outgoing particles. Akhiezer and 
Berestetskii present a convenient method to perform these operations. A 
brief outline of the method used follows . 

Assume that the cross section is given by da a |u 2 QuJ 2 , where 

u = u + y 4 . To find ckr a ^ ^ |u 2 Quxl 2 where ^ implies the 

Spl Sp2 Sp2 

summation over the final spin states of only one sign of the energy and ^ 

Spl 

is the sum over the initial spins of one sign of the energy, one first writes 


2 u (p) u (p) 

Sp 


m - i jar 
2E 


Now, 


I u 2 M u t 1 2 = u 2 M Uj ( u 2 M Ux) , 


since 


i u 9 u< i 2 = y u u y u * u 5 ' 
21 Li 2a la % 2/3 1/3 

a 8 H H 


- I 


a/3 


u U u * u * 

2a; la 1/3 2/3 


/v _i 


= u 2 U X u t u 2 


= u 2 Ux(u 2 Ux) + 




u 2 Ux = Uj u 2 = £ u la u 2a • 


a 


Thus, 


+ + 


|u 2 M u t | 2 = u 2 M \x t \i t M y 4 u 2 = u 2 M \i t u t M u 2 , 
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where 


M = y 4 M y 4 % 


One then notes that ( u u) . = u u 

'a/3 a 13 


u u = V u u =Sp uu 
LJ a a 

a 


Therefore, 


E E I u 2 Q u 4 1 2 

Spi Sp2 


= jr E E U 2 Q U 4 u x <5 u 2 

Spl Sp2 


= 4- E E s p [Q( u i u i) q(u 2 u 2 >] 

Spl Sp2 



where 


$ = Q(m - i^) y 4 Q + y 4 (m - i p) . 


Now, Q contains 0 = y e . To sum 
_ M M 

writes Q -0G, Q =0G, where G does not c< 

e . 

M 


I 
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Now fS = y 4 e y 4 = -& , since 

^ = e y + i e y 4 ; then = e y - i e y 4 and ee =e 2 -e 2 -e 2 = 1 
_ _ o'* ' o * /jl n o 

Thus, Q = -G0 and 

\ 2 |M fi |2 = "sFe" Sp G 0 (i& - m)} . 


Since definite polarization states of the photon are not of interest, 
one sums over the two independent directions of e^ . 

Let k be the space part of the photon propagation vector and the 
coordinate system be chosen such that the z-axis is along k. Thus, the two 

polarization vectors are taken to be e^= (1,0, 0,0) and e ^ = (0, 1, 0, 0) ; 

then, 

j- Z E Z ,M fi |2 = "iFE E SpL (j) G(ip' o -m) G^ (j) (^-m)l 

e Spl Sp2 o j=i ( f 


-We spjyjG(w- 0 - m > • 


(8-b) 


If in equation ( 8-b) , one of the polarization vectors e is replaced 

by kfi, thenX=y k = yik 4 + y 2 k 2 + y 3 k 3 + y 4 k 4 . However, k t = k 2 = 0 
M M ^ 

by choice, so k 3 = |k| , k 4 = i|k|; then y k = (y 3 + iy 4 )|k| , Therefore, 

Sp J(y3 + i 74) GU^-m) G#r(i^-m)| =0 

and 

Sp Jy 4 G(i^ o -m) G(y 3 + iy 4 ) (i^-m)J = 0 . 

If one multiplies the second equation above by -i and adds it to the 
first, one obtains 

S P jy 3 G(ip o ~na) Gy 3 (ijf-m) + y 4 G(i < p o ~m) Gy 4 (i^-m)J = 0 
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If, in addition, one multiplies this equation by -1/8E q E and adds it 

to equation ( 8-b) , one finds the summation carried out over the four values 
of j as follows: 

fZ Z Z i m £1 i 2 = -8Fe Sp K G(l V m) ■ 

e Spl Sp2 o * 


where v = i, 2, 3, 4. 

In this case, one desires to calculate the cross section 


, _ _e_ 4 1 V(q) I ‘ 
d<T 2k (2 tt) B 


7 E E I^Quil 2 


a e' 


-2- E E k 2 dk dfl. d£2 . (8-c) 


P o ° 


k p 


Therefore, the operation of summation over polarization states may 


be performed by the replacement I u 2 Q uj I 


F = |Q (ijar -m) Q (i^r-m)l , 
M o p 


8E E 
o 


SpF, where 


Q =Y 

jx r ix m i K i 


(if^-m) , (if'p-m) 

1 74 + 74 7, 


m/<2 M 


and 


+ (i-fi-m) 


( ii^-m) 


Q„ = 74 Q„ 74 = 74 y„ + y„ V 4 




m^! 'ju ‘jx m k 2 


Thus, 


Sp F = Sp 


(uf^-m) (i£-m) 

7 — ' 2 J 74 + 74 2 

r ix 4 m“K 2 


(i^ o - m) 


x 




(8-d) 


(ipr-m) 
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The spur of F in equation ( 8-d) can be written as 


SpF = F 1 (*c 1 , k 2 , q 2 , E q , E) + F 2 (/c lf k 2 , q 2 , E q , E) , 


where 


Fi =Sp 


*V if o" m) Y4 


(iii-m) 

mTq 


7^ 


and 

F 2 = Sp I Q (i^ -m) 7 y 4 (ija'-m) 

^ I M o M m 7 <: 2 4 

Since the substitutions P 0 ~*P> p-*P Q , and k-*--k in F 2 will 

result in an interchange of i<i-*k 2 , k 2 —K i , fi-»f 2 , and f 2 — *-fi, then 



F 2 (fCi, k 2 , q 2 , E o , E) = F 1 (k 2 , k u q 2 , E, E q ) . 


Therefore, one will calculate Fj only and then make the above interchanges 
to obtain F 2 . 


Making the substitutions for Q 


M’ 


one finds 


Ft = Sp 



(iij-m) 


74 +y 4 


( i-f^-m) 
m 2 K 2 


r, 


M 


x (i^f -m) 74 


(i£-m) 

mV 2 




+ 




74 (i P 0 -m ) 7 4 (i<[-m) 7^ (i^-m)j 

7 M (^ 0 - m ) 74( if r m ) 7^ ( i^-m) J . 
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The latter equation is a result of the fact that 


Sp(x + y) = Sp (x) +Sp(y) . (9-a) 

Other useful identities are: 

Sp CaOO = -J- Sp(jstff + J 2 te) =4ab ( 9-b) 


Sp (XXb'. . . ) = 0, if an odd number of operators is present ( 9-c) 
Sp [(Pi +m) (p$ + m)] = 4(pjp 2 + m 2 ) (9-d) 

Sp(xy) =Sp(yx) (9-e) 


Sp [AjA^sA^ - (AjA 2 ) (A3A4) + (AjA4) (A 2 A 3 ) 


“ ( A1A3) (A 2 A 4 ) 

y 4 b = -b* y 4 , b = any four vector 


y XMy = -2 jeW 


y ?C y = - 2gC 
M M 


y = 4ab 

V 


y^sy^ = 4s , s = any scalar 


7474 = 1 


y y =4 

V M 


77 7^ 7 = 4 6 „ 

' jj, ' ct r p r yL a(3 


(9-f) 

(9-g) 

( 9-h) 

( 9-i) 

(9-j) 

( 9-k) 

(9-1) 

( 9-m) 
( 9-n) 
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To calculate the spur of the matrices, first commute the y 4 matrices 
so that they are next to each other by using equation ( 9-g) . Next, use 
equation ( 9-1) to calculate y 4 . The sum over the gamma matrices is then 
accomplished through equations ( 9-k) through ( 9-n) . After performing the 
first operation, one has 


F 1 = ^2 (y^(iii-m) ( -i.gr* -m) (lij-m) y^(ig-m) 


+ m^7 2 ( Y4(i ^ 2_m) » 

where 

= ytf* . 

Next, one considers each spur separately. First, write F i“ F li + f 12 ; 
then 

Fl1 = ( y ju [if i 0 o ^ + mi i Fl ‘ + 2 im2r i 

- ^ - im 2 p* - m 3 ] y ■ ( ipr-m) j . 


Summing over the gamma matrices, using equation (9-c) to eliminate the 
spurs with an odd number of operators, and noting equation ( 9-a) , one writes 

F ll = ^T72 (sp [2Fi or* fipr] + Sp[4m 2 f i g] + Sp [m 2 fipi 

Hi Kj[ I O 

- Sp [ 2 m 2 ^* g] - Sp [4m 2 f 1 f I - 4m 2 f!P o * - 4m 2 p o *f 1 - 4m 4 ]| . 
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Equation ( 9-h) is now used to evaluate the spurs to obtain 

Fn= [zifiP* *iP) ' ( f i f i) (P Q * P) + 4m 2 fiP o * 

- 2ra 2 f 1 f 1 - m 2 p o * p + 2m 2 f 1 p + 2m 4 | . 


The next step in this calculation is to compute F 12 , where 
Fl2 = {74(^2“™) r^(i^ o -m) (i 74^1 -y 4 m) 7^ 


Following the same procedure as before, one finds 


F 12 = — — (-2-f^ ~ 2m 2 ff & - 2 m 2 i' 2 ‘' # 

12 m i K i K 2 \ o 1 ^o 4 * 

+ 4m 2 (p Q4 + f 14 ) (p 4 + f 24 ) -2m 4 )J 


in which has been added the relationship given by equation ( 9-n) . 

The spur can be computed by the previous method to obtain 

Fi2 I - (f * fi) (p 0 p) ■ (f 2 p) (f * p 0 ) + (f 2P c > ( f r 


P) 


- m 2 (ffp Q + f 2 "p) + 2m 2 (p Q4 + f 14 ) (p 4 + f 24 ) 


Therefore, for F 1? one has the expression, 


Fj = F u + F 12 = -f-g - { 2 f 1 P Q * fiP - fi fi P Q * P + 4m 2 f lP ^ -2m 2 f 1 f 1 

- m 2 p * p + 2 m 2 hp +2 m 4 } + — r- 
o 1 J m 4 K 


K\K 2 


X 


j - (f* fl‘ e ) (P Q P) - (f 2 * p) (f* P Q ) + (f 2 p 0 ) (ff p) 

-m 2 [ff P q + f* p + 2 (p Q4 + f I4 ) (p 4 + f 24 ) - m 2 ] J . 


( 10 ) 
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After performing the required algebra, one obtains 


T F * 

4 k { 


*1*2 


+ 


2 \ 4E 2 4E E 

2*) - (*i +2 > + -^ *i 


m 


m“ 


Kit 


q‘ 

~ 2 

+ 

q 

K2 '^ 

- 2 J + 

2 E 2 

U- 


2 E E 
0 

+ _ ^ r 

Tn 2 / 

m 2 


( 11 ) 


+ k 2 + 4) 


After interchanging Kj with k 2 and E^ with E one obtains F 2 ; 


i. e. , 


T F!= 7J | ( K * Kl + 1^) ' 


.„ 2 4E E 

4E 2 , o 

T («2 + 2 ) + — IX” k 2 


m 


m 


+ 


K t K 2 


•^(‘2 + «l-^ ' 2 ) + ^ 

« 

2 E 2 / 2 \ : 

o / q 2 \ 

m 2 \ 2 m 2 / 


F" («2 + + 4 ) 


I 

If — ^ Yj l"u 2 Q u-i 1 2 is replaced in the differential cross section, 

M M 

7 e 1 

equation ( 8 -b) , by ■ ■ ■■ £ (Fj+F^ , one has 

o 

fi E Epdfi d£2, 

, e° o dp dk , 2 p k 

I V( q) I 2 XT X 7^7 k 


4kE E E p (27 r) 

o o 


— 2 2 
«1 *2 


q 2 

“ 2 ^1^2 


2 

— 2 (E 2 +E 2 )+/c 1 + /Co 2 “ 2 

m o 1 2 m 




+ (Ki + Ki)[KiK 2 + ^J ~~2 (KjE + KgE^ 2 
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Now a = e 2 , and, as will be seen in the ensuing section, 

| V ( q) 1 2 = (Z 2 /q 4 ) [1 - F(q, Z)] 2 , where F(q, Z) is the atomic form factor. 

Since 

m 2 /c 1 = 2(p • k-Ek) = - 2k(E - p cos 9) , 


m 2 K 2 = -2({T • IT- E k) = 2k(E - p cos 9 ) 
* o o o o o 

q 2 = mVj + m 2 «2 - 2 m 2 + 2 E q E - 2 p^p , 

p 2 = E 2 - m 2 , and p 2 = E 2 - m 2 , 
o o 


( 12 ) 


one writes dcr as 


do 


Z 2 q 3 
(2tt) 2 


dk 
P o k 


dfi 



P 


[1 


F(q,Z)] 2 


x 


p 2 sin 2 9 
(E-p cos 0) 2 


(4E Q 2 -q 2 ) + 


p 2 sin 2 9 
o o 

(E -p cos 9 ) 2 
o o o 


(4E 2 -q 2 ) 


p 2 sin 2 9 + p 2 sin 2 9 

+ 2 

(E -p cos 9 ) (E-p cos 9) 
o o o 


(13) 


2 p p sin 9 sin 9 cos <b 
o o 

(E -p cos 9 ) (E-p cos 9) 
o o o 


( 4 E E - q 2 + 2 
o 



This is the Bethe-Heitler cross section for bremsstrahlung production, 
differential in photon energy and angles. The atomic form factor, F(q, Z) , 
is available only in numerical form; therefore, one cannot integrate equation 
(13) in closed form. However, from the discussion in the next section one 
can, for present purposes, set F(q, Z) = 0 . 

The momentum transferred to the nucleus, q, can be written in 
terms of the angles. Since q = p^ - p - k, then 
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q 2 = p 2 + p 2 + k 2 - 2 p k cos 9 + 2 pk cos 9 

o c o 


- 2p p (cos 9 cos 9 + sin 9 sin 9 cos d>) 
o o o 


( 14 ) 


The element of solid angle in equation ( 13) is 

dfi d£2. = sin 9 d 9 sin 0 d 0 dd> d <£ 
p k o o o 


The integration of equation (13) over the angles is straightforward 
but tedious. (The integration is indicated in Reference 16. The result of 
the integration over angles is [ 2 ] 

/ 

y « a \ 

eE 


, Z 9 dk P 

dcr = -tt- r -r 

137 o k p 


2 E E 

3 o 


P 2 + C 

7^ 


e E 

1 + ——rr + 


e e 


P o P 


+ L 


8 E E (e 2 E 2 + p 2 p 2 ) 
+ k? \ o o ; 


3P 0 P 


P o V 


(15) 


k 


2 p o P 


E E + p 
o o 

n — 


e ~ 
o 


E E + 
o 

— rs — 


2kE E 


where 


E = E + k , 
o 


L = 2 In 


E E + p p - 1 
o o 


e = 2 In ( E + p ) , 

o o o 


e = 2 In ( E + p) , and 5 - = r 2 

' v ' ’ me " 2 o 

o 
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Equation (15) gives the bremsstrahlung cross section differential in 
photon energy, k, and is the expression for which one wishes to obtain a 
Coulomb correction factor. Note that equation ( 1 5) contains k explicitly; 
thus, a spectrum of photon energies for a given incident electron kinetic 
energy is obtained. 

SCREENING BY THE ATOMIC ELECTRONS 

Equation ( 8-a) included a factor X( q) , where y\ q) is the Coulomb 
potential in the momentum representation. To determine the importance of 
including screening effects in one's calculations, this factor will be investigated 
further. 

One first writes the potential in the momentum representation by 
performing a Fourier transform on the interaction potential for bremsstrahlung 

[14]: 


y^q) = / V(r) exp(iq • r) dr , 

where r is the radius vector and q is the momentum transferred to the 

nucleus in the process. Note that V(r)=V +V , where V is the 

n e n 

potential resulting from the nuclear charge and is the potential resulting 
from the atomic electrons. 

If the atom is assured to be spherically symmetric, the total charge 
density can be defined as 


p(r) = p (rj - p (Vj , 

ii y 

where p^ is the nuclear charge density and p^ is the electron charge 
density with p^(r) > 0 and p^ (r) > 0 . Then, 

V(q) = ^7 / [p n (rj “ P e (rj] e 1Q r dr 

(16) 

[V*- F e<*] * 


31 


Now, F (q) « 
n 

1 , if p n (rj w Ze 6 (r) 

F e (q) = J 

r ; i q • r 

' p ( r) e dr 

e 

2tT 

Ze • 

r 2 j : i q r cos 6 

r drp (r e c 

J e 

_ 47T 

Ze • 

r , ? sin (qr) 

/ dr rp ( r) ‘ 

> e qr 

47 r 

Ze q 2 

/ M 2 P “ j (M) dM . 
J q o 


( 17 ) 


where p = qr and j (p) = Sm - . Thus, l^(q)l 2 is proportional to 

O fJi 

(Z/q 2 ) [1 - F ] 2 . If screening is not important for a given problem, then 

G 

F = 0 . 
e 


To determine the quantitative importance of screening in our energy 
and atomic number range, one first investigates the important impact 
parameters of the process. If a classical approach were applied to the 
problem, one would investigate whether the Coulomb field is screened to a 
great extent for the main contributing impact parameters [14] . Since a 
plane wave is involved, however, an exact definition cannot be given to an 
impact parameter. This is generally true in quantum mechanics. To give a 
rough indication, however, note that the integral in equation ( 16) implies an 
averaging over all impact parameters. 

If one examines the zero-order Bessel function in equation ( 17) 

( Fig. 4) , one finds that the main contribution to the integral comes from 
around qr < 1, since for qr > 1 the function follows a damped oscillation. 
However, for qr < 1, the factor r 2 in the volume element is small. Thus, 
the main contribution comes when qr « 1 or r= i/q. This is defined as 
b; i.e. , b = 1/q, the most important impact parameter for bremsstrahlung 
production. Now q will have a minimum value when the momentum of the 
electron, both before and after photon emission, is parallel to the emitted 
photon [12] . By equation (14) in the preceding section, it can be seen that 

lq . I = fp -p'-'kl = IPJ - iFl ' (E -E) . 

min o o o 
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If we assume E q » i, then E - Ip I « 1/2E or, upon temporarily restoring 
the units, 

Vin = (m o /2) < 1 /E - 1 /E o> 


» k (m c 2 ) 2 / 2E E . 
'o' o 



Figure 4. Amplitude of the radial function of the atomic form factor. 
Since b will be a maximum when q is a minimum, 

b = fi c/q . = (fi/m c) 2E E/(m ck) 

max Tnin o o o 

b « KE E/k . 
max o 

Thus, for a given k/E, b^ will be greater for higher incident electron 

energies. It would therefore be expected that for higher energies, the 
screening of the Coulomb field by the atomic electrons will be important. 
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Next the extreme case will be considered; i.e., b » a , where 

max 

a is the atomic radius. This is defined as complete screening [14] . If the 
atom is assumed to correspond to the Thomas -Fermi model, then 


-t/3 -1/3 

a » a Z = 137 Z Z ' 
o o 


where a is Bohr's radius of the hydrogen atom. Since b » a , one 
o max 

obtains the correct order of magnitude if one replaces by a, then 

l/3 

137 -X Z « 'i E E/k . 
o 00 


Therefore for incident electron energies, E^, on the order of 137 Z 
(k/E) , there will be complete screening. 


A comparison will now be made between the radius of the Thomas 

-i/3 

J-CJ I *-W I, 

TF 
If b 


Fermi atom (r, 
parameter 


max 


137 Z _/ in our units) and the maximum impact 
is much greater than the nuclear radius but less than 


the atomic radius and if one excludes screening effects, the differential 
cross section will not be overestimated significantly [ 5] . 


To investigate the specific conditions of interest here (i.e., incident 

electron energy from 0.1 to 2.0 MeV and atomic numbers 13 through 79) , 

one follows the technique used by Koch and Motz [ 5] and compares the 

maximum impact parameters and the Thomas -Fermi radius over a wide 

range of incident electron energies. (Figure 13 of the review article by Koch 

and Motz has an error in the ordinate units. The results are plotted in units 

of 10 cm instead of X as indicated. ) One first calculates b for 

o max 

these energies and for a spectrum of photon energies in units of the electron 

Compton wavelength, £ . The results of these calculations are shown in 

Figure 5 where r has been superimposed for several values of Z. It 
Tr 

should be noted that for lower energy photons, screening is important (electron 
deflections occurring at a large distance from the atom) . For incident 
electron energies in the range 0. 1 < T q <2.0 MeV, screening is not too 

important except for low energy photons in high Z materials. However, as 
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can be seen in Section HI, the Born approximation is relatively accurate at 
the low frequency end of the spectrum when compared with experimental 
results. 

Therefore, since the primary interest here is Coulomb corrections to 
the Bethe-Heitler theory and it has been demonstrated that in this area of 
interest screening is not too significant, it will be assumed that ^(q, Z) = 0 

and that the derived correction factor will correct for Coulomb effects of'the 
Born approximation only. 

D. The Bremsstrahlung Problem Using Coulomb 
Wave Functions 

In Section II. C the electron initial and final state wave functions were 
assumed to be plane waves. This resulted in the Born approximation for 
bremsstrahlung cross sections. The theoretical results, at moderate 
energies, may be improved significantly by using Coulomb wave functions 
instead of plane waves in the matrix element. 

The second formulation of the problem will be approached by first 
separating the Dirac equation into polar coordinates (Appendix C) . The 
radial equations will then be solved for a pure Coulomb potential. The matrix 
elements will next be obtained by partial wave expansions for the incident 
and scattered electron and for the photon. Finally, the expression for the 
differential cross section will be derived in a form applicable to a computer 
solution. 

1. THE MATRIX ELEMENT FOR BREMSSTRAHLUNG PRODUCTION 

The Hamiltonian of the Dirac equation in the case of electromagnetic 
interactions, as is the case for bremsstrahlung production, is 

U. = c~a_ • ( fT - X) + §_ + V , 


where A is the vector potential of the radiation field and p, a _ , and £ 
are defined in Appendix C. 

The term - e a • A in the Hamiltonian is the perturbation responsible 
for the transition of the system in the bremsstrahlung process. 
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The matrix element for bremsstrahlung production is then 


M fi = _e f r r’(*f + “ ^) dS 


( 18 ) 


where # (^) is the initial (final) electron wave function. The final state 

wave function may be written (Appendix C ) : 




r S K' 


ir 1 f £2 


(Throughout this and the following section, the unprimed variables and 
quantum numbers refer to the incident electron, while primed variables 
refer to the scattered electron. ) 


The initial electron wave function has the asymptotic form (r 
of the superposition of a plane wave and an outgoing spherical wave 


00) 


ill ~ 

i 


u(m) e 


i k 


r , -l ikr 
+ a r e 


where a is the amplitude of the scattered wave and m is the z component 
of spin associated with the plane wave. 


To obtain a calculable form for the initial wave function, one expands 
the wave function into spherical waves. The expansion of ip in spherical 

waves, normalized in the energy scale, can be shown [17] to have the form 


Ip, = 47T 

i 



Z 


1+1 

Z 

m 


l 2 j 

s c K K 
k H - m, m 


, r * M - m r , M 

x Y / (P) ^ 
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where 



and s are constants, the values of which are determined from the 

K 

asymptotic form of the radial functions. In the Coulomb case, one finds that 
the difference between the Coulomb phase shift and that of the plane wave 

i 6 ' 


is given by 



m 


X¥ i 


s = e 

K 

K 

(Appendix B) . For ip-, may write 

(* \ 

i 1/2 

id' l +1 S. 5 j 
Y Y e K i K C K K 

\ 2 E p y 

! 

u u u - m, m 

KH m 

c fji - m 

(P) 



This can be simplified by taking the direction of propagation of the incident 
electron along the axis of quantization; say, the z axis. One first writes 
the spherical harmonic as 




M 


(cos 0 ) 
o 


im q> 


M = jit - m , 


and from the recurrence relations for 



one has [ 18] 


pM + i 
i 


(x) + 


2M 

(1-x 2 ) 1 / 2 


P_ M (x) + (i+M) (i-M+i) 

Xj 


p M-l 
i 


(x) =0 


where x = cos d . For propagation along the axis of quantization, 

0 = 0, x = 1, and the recurrence relation becomes 

o 
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p M+l 
i 


( 1 ) + 


2M 

(i-i ) 1 / 2 


P i M (i) + (£+M) (i -M+l) P^" 1 (i) 


0 . 


The first and last terms are finite but the second term becomes infinite as 

M M 

x — i unless M = 0, P^ (1) = 0, or both. That is, either = 0 or 

M = 0 . One must choose M = 0 (p = m) ; otherwise, the wave function would 


vanish. Now [19], Y^° 

one writes Y^° (0 ,cp) = 



i/2 

P^ (cos0) ; and since P^ ( 1) =1 , 


After substitution and simplification, this expansion becomes 

ip = — 

m ( 2 E p ) 1 / 2 
o o 


E E 


m 


iS'f+i i gj 

e * i K (21 +1) 1 / 2 C K K 
k 0, m 


x 



Expansion of the Electromagnetic Wave 


The electromagnetic wave denoted by the vector potential of the 

radiation field, Aj— in equation ( 18) , can be written as a linear combination 
K 

of waves that are circularly polarized [19] . A One first assumes that the wave 
is propagating along the z axis. Thus TT=k |k| and one may write 


A r (z) = X a A ( 19-a) 

k P=±l 

where, by definition, ap = + 1 for a left-hand circularly polarized photon and 
ap = - 1 for a right-hand circularly polarized photon. 
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Following the method used by Rose [ 19] , the electromagnetic wave 
denoted by in equation (19-a) is expanded as follows. For a plane wave 

propagating along the vector k one writes the vector potential as 

— f- A 1 k * T* 

A=Nee , (19-b) 


where N is a normalization factor to be determined and e is the polarization 
vector for the photon; i.e. , a unit vector perpendicular to the propagation 
vector k and pointing in the direction of the electric field vector. 

One then writes 

A 1 A . ^ A v 

e = (e 1 + iPe 2 ) 

n/T 


with P = +1 (-1) for right (left) circular polarization by definition. 

If the z axis is chosen to correspond to k, then = e^, e 2 = e , 

e 0 = e , where e , e , and e are the unit vectors along the x, y, and 
z X J z 

z axes . The vector e may then be written in terms of the spherical basis 
vectors [19] where the vector B in the Cartesian basis can be transformed 
to the spherical basis, B’, by B' = U B, where 


U = 



or 


B , (B + i B ) 

+ nTT X y 


B 0 = B 


B = 


\ J ~~2 


(B - i B ) 

x y 
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Therefore one writes e in terms of the spherical basis 


A 

e = 


-P 





where 


are the spherical basis vectors 


„ _ i a . v 

k A - + (e ± 1 e ) 

nTT x y 


A 



A plane wave may be expanded by using the well known Rayleigh 
expansion; for an example, see Reference 19. 


i k • r = 4 tt 7 i in a ‘ ' m a 

3 « L i j„ (kr) Y (r 2 ) Y„ (r«) 

i , m v ’ l u l v v 


where j (kr) are the spherical Bessel functions that are defined in terms 
of the regular Bessel functions as 

j < (kr) “(zFr) J l+i Ur) • 

From the spherical harmonic addition theorem, 

V m -A * m ,a . (21 + i\ V 2 _ . 

2 Y , < r >> Y , (r ’> = Hsr) V oose) • 

m x ' 


one may express the plane wave as 



However, 

has 


(cos 6) = 




A 


(r) 


then, after substitution, one 


ik z 
e 


(4tt) 1 /2 £ (_i) £ (2i+l) 1 / 2 j (kr) Y 0 (r) 


1=0 


( 20 ) 


Therefore, 

A p = -P ? p N ^ ("i/ [47r(2i+l)] 1 / 2 Yi ° (r) j (kr) 
i 


The normalization factor, N, will next be determined to complete the 
expansion. 


One first writes the vector potential representing a classical plane 
wave in terms of time and position; i. e. , 


A("x, t) = N n e 


i(k • x - (a t) 


A 


where N is the normalization factor and n is the direction of propagation. 
This may be written as 


A(x ’ t) = A REAL (x ' t)+A IM (x>t) 


= N n [cos (F • x - w t) + i sin (k * x-cot)j 


= N n [cos (F * x - co t)] 


Following the method presented by Feynman [12] , A is 

xi 

normalized to give unit probability per cubic centimeter of finding the photon. 
The average energy density is tf.cc . 
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From Maxwell’s equations, the electric field is 


t 1 BA — 

E = - — - -- - V <£ ; or, since 0 = 0 , 


E = - 


c 3t 


Nnco . — _ 

S1Q ( k • x - co t) 


Now |B| = |E| for a plane wave. Thus, one may write the average energy- 
density as 


< s >av - i < |E|!+ ir ' ! >Av * i < IE|! > 


AV 


or 


■\ 


2 , , 2 


N co 
AV 47r c 


j sin 2 (Ic • x - cot)^)> 


N 2 co 2 
AV ~ 87TC 2 


since 


< 

We> Av - -j ■ 


However, 

units 

■\£^> av ; then, N=( 

8?rc/ “ J , and in dimensionless 

^ CO / 

T— ( 

li 

o 

1! 

, N = ^‘ /! . 

CO 


Thus, 




REAL 


/ 87t\ V 2 a 

l — J n [ cos ( k • x - co t) ] 


or 


V(x,t) = (~) / n ■— | exp [ -i ( k • x - cot) ] + exp [ i ( k • x - cot) ] 
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By taking the minus exponent, since one is considering the emission 
of a photon, and noting that a) ~ k , in our system of units one has 



where e is the polarization vector. (See equations (19-a) and (19-b) .) 

Finally, using equation (20) and the above normalization factor, the 
expansion for an electromagnetic wave propagating along the z axis is 

Ap=- p f p (x) l/2 Z (V [47r(2i+i>] 1 / 2 Yi °(r) ^ (kr) . 


Next one defines the "vectorial spherical harmonic" as [19] 


T s y c iix Y m '~ v t 

Aim 1 t-t m T -v, v i S 
v 


However, 


j ^m 2 m^m j t J, 


has the inverse expansion 


( 21 ) 


m 2 = y c jij 2 j ^ 

Ji h Y m 1 m 2 m r j 


This implies that it is possible to write equation (21) as 


2 

A 


,1 1A 
'OP 


Ai P 


where A is the parity operator, A = i , i ± 1 
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The expanded vector potential may now be written 


Ap= “P (y) 1/2 (4tt ) i/2 £ ( -l) 1 ( 2i +1) 1 / 2 (kr) c“ A T 

J? , A 


Xi P 


To obtain a photon propagating along an arbitrary direction, one 
utilizes the rotation operator [19] , which has the property 


RT 


(k) _ 


+k 


- 2 


q'=-k 


(k) (k) 

q'q q’ 


( 22 ) 


The vector potential for a plane electromagnetic wave propagating 
along an arbitrary direction may then be written 


r— 00 1 + 1 +\ 

= 2 2 (V (2i + D‘A 


n/F f =0 A=X-1 m'=-A 


x i (kr) C*** D X ,_(a>, 0, 0) T „ , 

J i OP m ! P r k k Xi m T 


(23) 


From equation (18), it is noted that 


M = JVL. 
fi 




a ib ) 

= m 


The matrix element M(P) , 


since 



- z 

p=±l 


a p A p 


where M= ^ a M(P) , 
P=±i 


will be considered 


Since the purpose of this section is to derive the bremsstrahlung 
cross section differential in photon energy, the angular or polarization 
details of the scattered electron are not of interest. Therefore, one may 
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express the final electron wave function as 




r "‘ g «’ a KV 


ir *f , S2 . . 
K -fC'fX' 


If an investigation of the angular distribution of the scattered electron 
were within the scope of this paper, one would expand in spherical waves. 

These waves would behave asymptotically as plane waves plus convergent 
spherical waves [11] with an expansion of the form 


. i +1 -i 6 t 

ip (E) = 4tt 3 / 2 (2Ep) ^ Yi i K e 


m 


K ! jw T m M 


r g K' 

xcf 7 m „ <^"( 8 ) ' 

Ai'-m",m" i . _ 

' ir f Q 

k' -k'ju 


2 ^’ .^ju’-m" a 


Now one substitutes equation (23) into equation (18) and the matrix 
element becomes, after equation (21) is substituted for T p : 


M(P, m . K ’. ft ') = - e J d^x 


- P 2 ^ ZEE <-»' 


°0 jg+i +\ 


.1 


n/T i=0 A=i-i m T =A 


X j, < kc > <f k V» Z l 


(24) 




iir f , , , 

K 1 -K'M 


2 2tt (2E o P o )" l/2 i K (2^+i) 1 / 2 

K 


r g Q 

, | . . * , b K KH 

,s , 

X Om 6 \ . -1, 0 

ir f fl 

K -KH) 
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Equation (24) will next be integrated over the angles of the volume 

element. The square of the absolute value of the result will then be put in 

the form |M„ I 2 = M* M„. . This expression will then be used to obtain 
fi fi fi 

the desired differential cross section. 

The integral in equation (24) may be written as 


1(1, m V. k.k'.M.M') = /d 3 xj (kr)Y m V (r) 


(25) 




lr *f (a ) £2 

K-V “KM 


r g (<7 ) fi 
fC = ^ kM 


A _/|\ A _/0 l\ / a 
since tf L) = VMl 0 


b r S Aa 


where j| and b are any two matrices. Thus, 


r Q m T - v A 

1(1, m ',v,k,k',h,h')= J d®xj (kr) Y (r) 


—2 + 

ir g . f fi , , cr fi ,, 
& k' k -k'm =v -km 


- 2 . 


i r f , g fi , , <7 fi 
k' & k -k'm' = v km 


s Ii “ I 2 


One also defines the radial integrals as 


K(k'k) = 

j[ (kr) g K ' 

f dr 

K 

K(kk') = 

OO 

X 4 (kr) f K' 

g dr 

K 
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then. 


I, = KU’k) / do (i) Y i m '- V <?)[ 0^, g„ 0. w ] 


(26) 


The quantity in square brackets in equation ( 26) may be expanded as 


follows : 


i 1 i 

r + ^ 1 v K ' 2 V + 

fl . . <7 Q = ) , C , Y 77 <7 

l k'ix'-v -khj u /x'-m,m i f m = 1 


a 5 j 

£ c -k^-k /*-« 

Lj u-a , a SL a 

a -k 


Now, 


+ 

Tj (J T} 

m -v a 


= (<j ) = 2 /j m | J I i a\ 

V p/ma \ v / 


From Reference 19, one has the expression 


then, 


(%) 


= 2( -1) V 6 C 2 1 2 

ma m, a+v 2 a+v,-v 


Therefore , 

[V, , a 0 1=2 (-l) m -“^3cV iV 

L k'ijl' =v -kmJ _ /x'-m,m 

* pt'-m 


m, a; 


xC K K C 2±2 Y 
\x-a,a. m,m+a; f t 

K 


r li-a 


-K 


(27) 
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Now the integral over the angles can be evaluated by using: 


* m 3 v Y 

h h * 


/ dfi y" Ui3 y/“ 2 Y. 




(21 1 +l)(21 g +l)l */ 2 W 3 c £ 1 i 2 i 3 
4ir ( 2i 3 +l) 0 0 ’ 


then, 
/dS2 Y 


*/i’- m y^ _cv Y m,-m+Q! 


-k 


( 2i+i) (2! _ +1) 

K 

47 r ( 2i +1) 
k' 


1/2 


(28) 


Thus equation (26) becomes 


ii a , i , 

x C * C “ K K 

0 0 m T -m+cy , fJL-a 


I 1 =K(k'k) Yj ( _1 ) 
m, a 


m -a 


i f i i i \ j 

x c K K c~ K ~ K 

m \i-oi , a 



(29) 


Note that the only difference between the angular portion of Ij and I 2 
is the changing of k' to -k* and -k to k ; then, 


I 2 ~ K (k «') Y (i) (- 1 ) 1 

j8 , a 


[X 

~ ( 2i+i) (2i +1)" 

K 

/ 47T 



1/2 


x C 


_ K » 2 J_y.| 


i. J J. 


(30) 


K' g K K 


X 1 1 

C 2 2 


ii i , 

„ K -K 


ii i 


K -K 




ju’-m, m n-a , ol m,m+Q? 0 0 
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Now, define 


A( A, f m', a , m) = ( -i) (2i+i) £ ^ 

m, Q! ,K 

l+l £ b j i 6 

xi K {21 +1) 1 / 2 C* * e * [Ij -I 2 ] 


/~3 . _ . , V i A . K 

A = n/^t ^ 2i+1 ) ^ C m , -m+Q! , m-a 1 (2 £ +1) 1 / 2 


m, OL , K 


£ 5 j i 6 , 2 •> t / +1 

xC 0m 6 (1) (_1) K(kk) V-m.ml iF^T 


i 2 j ! , , £££,£££, 

xC ’ K C 5 ^ C ~ K K C ~ K K 

ix-a,u m,m+a 0 0 m'-m+a ,n~a 


- K(kk') 


2i +1 N 1 / 2 £ ij , i |j 

K \ „ ~K -/C T „ K K 


2i +1 
-k' 


[i 1 -m , m n-a,a 


Hi f ii i f 

c * "* c * ■* 


m, m+a 0 0 m'-m+a ,/i-q: 


Therefore, the matrix element, equation (24) , can be written as 

00 i +1 _ r-— 

E v* V V 2tt e \l 3ir 

L L L 

£ =0 A=i -1 K a , m \/k E p 

o o 


xPC op A D^,p (<P k 0 k O) A(AiK'M’ma) 
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and 


|M(P, m./c'M') I 2 =M*(P , 1 m ) K',icM')M(P,m,)c'K/J') 


1 2tt 3 e 2 

kE p . 

o o Xji i 


y P» D ^ ( (p 0 0) 

Lj %P» p’-m,P' v ^k k ' 


(32-b) 


xA*(X 1 ,! 1 ,K , ,M , ,m) X Z p C o 2 p X2D ^- m> p <<P k \°) 

^2 h 


X A(\ 2 ,h>«' >1*' m ) • 

These matrix elements will now be used to find an expression for the 
differential cross section for bremsstrahlung production. 

2. THE DIFFERENTIAL CROSS SECTION 

Recall from equation (7-b) that the differential cross section is given 
by 


2 ?^P E 0 E ^ IM fl l 2 ^dkd! 2p dS2 k . (33-a) 


M A* 

y e 


The density of final states (equation 6) included in equation (33-a) is 


©' 


p f = 4 (^) E p k . 


This can be written as the product of the density of final electron states and 
final photon states; i.e. , 


P f P f P f 
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However, one assumes, for present purposes, that the photon is 
observed and the scattered electron is not observed. Therefore, since the 
Coulomb wave functions have been normalized in the energy scale, one 
writes the density of final states as (since p f = 1 ) : 

e' 


P 


f 


= 2 



k 2 


and equation (33-a) becomes 


dcr " t 2 t lM fi |! k2dkdn i 

p ’ v o 

y 


or 


, E 

dCT ' = Z 7TT2 — I M | 2 k 2 dk dS2 

LJ p fi 

H O 

y 


k ‘ 


(33-b) 


The matrix elements, M(P, m,K',ju') , included in equation ( 32) are 
related to the matrix elements of equation (18) by the expression 

M fi = Z a p M(P,m,«: , ,M') • 

iL P=±l 


Therefore, the differential cross section is given by [11] 


<b '' E ( ,a P |!d<, i 

P=±l l 


PP + a -P a P d<T -PP 


(34) 


where 


E 

dCTp,p=(27r) 2 — Yj M + (P'mK'p') M(PmK'/Lt') k 2 dk d(cos 0 q ) d cp ° . 

P o mK'M’ ° ° 35) 
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Using equation ( 32) for M (P'm/cV) and since the 

azimuthal angle of the photon is not of interest here, one writes 


^prp = f k dk d(cos 6 ) d 0 


? P ' C 0P- X ' V-m,p.^o 9 o°> 


Mi 


x A' fi 


(36) 


y p c h 1 a 2 n A 2 

^ C flP r» 

M 2 


OP D M'-m,p (<p o ,0 o O) A (M 2 «Vm) 


To perform the integration over the angle, the integration over the 
rotation operators is first considered; i.e. , 


I = f d vAi 
■J p'- m,p* 


(<p e o) 

^ o o 


D 


/i'-m, P 


(<p e o) 

^ o o 


d cp 


Since the rotation operators are independent of <p , one may immediately 
write . 


I = 2tt D 


H'- m, P' 


A 2 

'-m, 


P 


However, 


M'~ m,P' 1 ; -M'+m, -P' 


Then, 


I = 2tt (-1) M ’ m_P ' D Al , D As 

-M’+m, -P» ju'-m, P 
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Now, since 


’'•i n *2 _ y P^-i^j n i 

-ju'+m, -P' M'-m, P L -At'+m,/x’-m -P’,P 0, -P'+P 


and -ju'+m+ju'-m = 0 , the integral becomes 


I = (2vr) (-1) 


M'-m-P’^ c Aj A 2 j 


•1^2 3 p^l ^2 J tJ 

p'+m,n'-m -P\ P 0, -P'+P 


Furthermore [19] D 


= 

\2i +1 / 


47r \ 1//2 * m 

£+1/ Y i (0,<P) ; 


then, 


I = (2tt) (-1) 


M'-m y p^l^j p^iA 2 j ( 4tt \ 1 / 2 P-P* . 

^ C -^+m,(U'-m C P'P (,2^+1 j Y j * 


Since [19] yJ m (0,«p) = (-1)™ m (0,<p) , 


one writes 


* p-p’ P'-P 

Y (9,<p)=Y (B,<p) . 


P'-P 

(The fact that (-1) = 1 has been used.) 


The spherical harmonic can be written as 


yP ,_p r2i±L] r ;j-;p;-p;;; l 1/z P i p, - p ' (008 9 , , 

J 47 r (]+|P'-P|)'. ] o’ 
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where 


p IP'-Pl (oose )a _J_ 

> ° 2 J j! 


(i~COS 2 0 ) m / 2 


d(cos 0 ) 
o 


j+m 


(cos 2 0 -i) 1 
o 


are the associated Legendre functions. Thus, 


1 = 2* 2 C Xl *l ] , C XlX2j 

L-t 7 -£i T +m .u’-r 


(i- ip'-pI) : 


(j+ Ip’-pI) i 


-/Li'+m,/n'-m -P',P 


V 2 IP' -Pi 

P. 1 (cos 0 ) 

3 o' 


(39) 


The differential cross section, equation ( 36) , may then be written 
(since a = e 2 in nuclear dimensionless units) 


2 

da^.^= — 2 ° kdkd(cos0 ) Y 


P'P p 


3 


(i- IP'-PD: 
(j+ IP'-PI)! 


V 2 | p* -p | 

P 1 (COS 0 ) 

3 o' 


V pp, V pM^2 P ?4*23 y p?4^2 3 

,1 k C 0P °P p,p , -M'+ m./u'-m 

X l Xs M (40) 


*1*2 


x ^ A* (X^jK'jU'ni) A(A 2 i 2 x'lJi'm) 
k ' m 


For unpolarized photons, equation (34) reduces to do ' = 2dcr{ > t . 

(The factor 2 has been included in the density of final states. ) Using 
equation (40) for doj j , one obtains the cross section for unpolarized 
photons 
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dcr ' = 


6tt 2 a 


rr- dk d ( cos 


x C 


■* 1^1 r .^2 * ^2 


01 


01 





C ^1 ^2 j 

-M'+m.jtt'-m 


( 41 ) 


x Yj A* (XjijK'jU'm) A(\ 2 ^2 K 'l jl ' m ) 
k ' m 

This may be simplified by considering the orthogonality relation for the 
C-coefficients 


y c ^i^2j c ^i *2 i _ « 5 

V -jit’+m, /P-m 1,1 i,-jU'+m i,ju'-m 


Therefore, equation (41) becomes 


d<7 ' = 6 7T 2 C2 



dk d(cos 9 ) Y 

j 


P. (cos 9 ) 
J o' 


x 2 

•M2 


2 

Ai A 2 


C 


illAi 

0,1 



(42) 


The purpose of this report is to obtain the bremsstrahlung cross sec- 
tion differential in photon energy. Therefore, one integrates over the photon 
angle, 9 . Note that the only angular dependence comes from the 

associated Legendre function; one thus considers the integral 
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Now, P (cos e > = i ; therefore, 
o o' 

I = f P. (cos e ) P (cos 6 ) d (cos 6 ) 

J j o o o o 

But, from the orthogonality relations for these functions, 



2 

2j + 1 


6 . =2 

J.o 


The only value of j that contributes to the summation is, therefore, j = 0 . 

One then has the bremsstrahlung cross section differential in photon 
energy, for unpolarized incident electrons, and without regard to the 
scattered electrons or polarization of the emitted photon 



Z E 

ilf 2 A. 1 x 2 


-,^1 1 r .^2 f ^2 


Y (-1 ) M Y A *(M 1 K , ,M , .m) 


Oi 


oi 


H' 


k 'm 


x A(A 1 i 1 «: , jU , m) 


(43) 


where A* (A 1 ,£ 1 ,(c', M'm) A( X 2 ,f 2 , /c’.jU'm) is given by equation (31) . 


SECTION III. CROSS SECTION RESULTS 


A computer program was developed by the author to calculate the 
Bethe-Heitler cross sections for any desired incident electron kinetic energy, 
a complete spectrum of photon energies, and the scattering media of interest. 
Results were initially obtained for incident energies from 0. 1 to 2. 0 MeV 
for aluminum, copper, tin, and gold. These materials were chosen so that 
available experimental data could be utilized for comparison purposes and 
so that a wide range of atomic numbers would be used in the calculations. 

The results of these calculations are shown in Figures 6 through 9. 
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Figure 7. Bremsstrahlung cross sections differential in photon energy for 0. 5 
MeV electrons incident on aluminum, copper, tin, and gold. 








The experimental data for 0.2, 1.0, 1.7, and 2. 5 MeV were obtained 
from Dance, et al. [7] while the data for 0. 5 MeV were obtained from Motz 
[6] . It is noted that the experimental results of Motz give, in general, a 
higher value for the differential cross sections than the results of Dance, 
et al. , who attribute this difference between the two experiments to the 
presence of background in the experiment of the former. Since experimental 
data between 0.2 and 1.0 MeV were required for comparison purposes, the 
author had to use the 0. 5 MeV data from Motz as the only available data in 
this range. 

The computer program first written by Zerby and Brysk [11] and 
expanded, reprogrammed, and improved by the Space Sciences Laboratory 
of NASA's Marshall Space Flight Center was used to calculate the 
bremsstrahlung cross sections using Coulomb wave functions. These results 
are also included for gold and aluminum. 

It is noted from Figure 6 that for low incident electron kinetic energies 
the deviation of the Bethe-Heitler theory from the experimental results is 
appreciable. This deviation increases markedly for higher atomic numbers, 
as expected. In fact, it can be noted from Figure 6 that the theoretical 
(Bethe-Heitler) cross section for gold is actually below the experimental 
results for tin as the high frequency limit is approached. However, the Bethe- 
Heitler results improve when compared with experiment as the incident electron 
kinetic energy increases, as expected. 

The broken lines in Figures 6 through 8 indicate the differential cross 
sections that were calculated by using the correct Coulomb wave functions 
and the partial wave expansions. Note the very good agrement, in this energy 
range, between the experimental results of Dance, et al. and the theoretical 
calculations. 

Table i illustrates the importance of including the Coulomb effects 
for the energies that are considered here- Table i lists the relative error of 
the Bethe-Heitler theory because of the exclusion of Coulomb effects, assuming 
that the partial wave expansion approach results in the correct differential 
cross section. 

From the discussion in Section II. C, subsection entitled Screening 
Effects by the Atomic Electrons, one can see that the small deviation of the 
Bethe-Heitler theory above the experimental values at the low frequency end 
of the photon spectrum (Figs. 6 through 8) is a result of the exclusion of 
screening in making the calculations. In addition, note that the differential 
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TABLE i. THE APPROXIMATE RELATIVE ERROR OF THE 
BETHE-HEITLER THEORY WHEN COMPARED TO THE 
PARTIAL WAVE RESULTS 



\ k/T 







tN. ° 

o \ 

0.2 

0.4 

0.6 

0. 8 

0. 9 


0.2 

0. 05 

0.26 

0.48 

0.65 

0.76 

Z=79 

0.5 

— 

0.23 

0.41 

0.61 

0.75 


1.0 

— 

0. 18 

0.37 

0.53 

0.62 

: 








0.2 

0. 02 

0.04 

0. 16 

0.29 

0. 54 

Z=13 

0.5 

— 

0. 01 

0.12 

0. 16 

0.43 


Relative Error 



cross sections calculated by using the unscreened partial wave approach are 
also too large at the low frequency end of the spectrum. The much larger 
deviations of the Bethe-Heitler theory at the high end of the spectrum are 
almost entirely a result of Coulomb effects. This is shown schematically in 
Figure 10 where the differential cross section is plotted as a function of the 
ratio of the photon energy (k) to the incident electron kinetic energy (T ) . 


SECTION IV. COULOMB CORRECTION FACTOR 
A. Extension of the Elwert Factor 

The general approach used to obtain a correction to the Bethe-Heitler 
theory was to first calculate the cross sections, differential with respect to 
photon energy, using both the Born approximation and the partial wave 
expansion approach. Next, it was assumed that 
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k/T 

o 


Figure 10. Schematic representation of screening and Coulomb regions. 



= f(lM 0 .Z) 



(44) 


where 



is the differential cross section using the partial wave 
expansions; 



is the differential cross section obtained from the Born 
approximation; 
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0^ = Pq/Eq » /3 = p/E as previously defined; 


Z is the atomic number of the scattering medium; 

and 

f(/3,/? o , Z) is the correction factor to be derived. 

That is, it was initially assumed that the correction factor would be a function 
of the speed of both the incident and scattered electrons and the atomic number 
of the scattering material. The correction factor may then be written 


f(/3,0 o ,Z) = 



(45) 


To get an indication of the analytic form of this factor, the ratio of the 
cross sections was first computed and plotted as a function of the ratio of the 
photon energy (k) and the incident electron kinetic energy (T ) , as shown in 

Figure 11. From the resulting curve, it was first noted that the factor was 
similar to the Elwert factor in form (Fig. A-l(b) of Appendix A) . There- 
fore, as an initial approach, it was assumed that the factor has the form 


0 Q { 1-exp [-g(Z)/0 ]} 
f(0,0 o ’ Z) = 0 {i-exp[-g(Z)/0]} 


(46) 


and determined the function, g(Z) . 

Using multiple regression techniques, based upon a least squares 
criterion, values of g(Z) were obtained using gold and aluminum as the 
scattering media. The most accurate function, assuming that the form given 
by equation (46) is correct, is shown in Figure 12. The solid curves are the 
actual correction factors obtained from equation ( 45) while the broken line 
curves are from the regression equations. It was noted that the results are 
not very accurate; therefore, this analytic form for the factor was abandoned. 
It is of interest to note in passing, however, that, at an incident electron 
energy of 0. 2 MeV, the results of the regression equation are almost 
identical to those of the Elwert factor. 



5.0 


Z = 79 

Z = 13 



0.2 


0.4 


0.6 


0.8 


1.0 


k/T, 


Figure il. The ratio of the differential cross sections calculated by the 
partial wave method to the Bethe-Heitler theory. 


B. Alternative Forms 


From Figure 11, one can see that the curve has a value of 1. 0 at 
k/T^ = 0 and has the general shape of a hyperbolic cosine function. Since 


a cosh (x/a) = a/2 [exp(x/a) + exp( -x/a)] , 
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ACTUAL 

COMPUTED 



Figure 12. Extension of the Elwert factor to higher energies for incident 
electron energies of 0.2 and 0.5 MeV incident on gold. 

it was assumed that the correction has the analytic form 


f 


2 


J exp |Vg (j 8 q , Z)J + exp £-y/g(, 8 q , Z)J 


(47) 


where y = k/T^ and g(P Q > Z) 


is an unknown function to be determined. 


6 




A curve fit for f in terms of g(/3 Q , Z) was first obtained. This 

resulted in different values of g for a given atomic number and specified 
value of /3 q . The results are shown in Figure 13, where the regression 

curve is compared with the actual curve. 



Figure 13. Comparison of calculated to actual correction factor. 

The second step was to obtain the analytic form of g(P 0 > Z) . It 

can be seen from Figures 6 through 9 that the magnitude of the correction 
should increase with Z and decrease with increasing /3 . Therefore g 

was plotted versus the parameter a = jS^/Z , since the exponents in 
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equation (47) are defined as y/g . Using standard regression techniques 
[20] , the following function was obtained: 

g( j 8 Q ,Z) = 0.394+ 9.47 (j8 /Z) . 


Slightly more accurate functions were obtained but these had more 
complicated forms. Therefore, since the linear relationship gave relatively 
accurate results, it was used in these calculations. 


The final form of our correction factor was found to be 




" k 


k / i \ 


2 

exp 

_ T o ' g /_ 

+ exp 

_' T o \Z> _ 



( 48-a) 


where 


g = 0.394 + 9.47 (j3 /Z) 


This can be written 


f = cosh 


— ( - 
T I 0. 

o \ 


394 Z + 9. 47 (3 


:) 


( 48-b) 


C. Comparisons With Experimental Cross Sections 

The corrected cross sections were next calculated, using equation 
(48-a) in conjunction with equation (44) , to make a comparison with the 
experimental cross sections. The results of these calculations are shown in 
Figures 14 through 16. 

An inspection of Figures 14 through 16 indicates that the corrected 
Born approximation gives relatively good results for the energy range of 
primary interest. However, it may be noted from the 2. 5-MeV curve in 
Figures 14 and 15 that the correction factor does not accurately predict the 
transition of the Born approximation from below to above the experimental 
results. This occurs between 1. 7 and 2. 5 MeV and was discussed in 
Section I. The failure is particularly apparent for high-Z materials. It is 
therefore concluded that the uncorrected Born approximation gives more 
accurate theoretical results in the energy range above, say, 2. 0 MeV. 
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(cmVMeV) 









It may be noted from Figures 14 through 16 that the corrected 
differential cross sections are closer to the experimental results for lower 
energies. This is most likely because of the fact that the partial wave data 
used in the derivation of the correction factor were limited to lower energies. 


SECTION V. CONCLUSIONS 


The bremsstrahlung cross section has been calculated using both a 
partial wave expansion that includes correct Coulomb wave functions and the 
Born approximation. It was found that the former approach gives better 
agreement with experiment in the 0. 1- to 1. 0-MeV energy range. In fact, 
the Born approximation differs appreciably from the experimental data, 
especially for high atomic numbers and lower energies. For example, the 
greatest deviation was found to be about 75 percent for 0. 2-MeV electrons 
incident on gold if the photon received 90 percent of the incident electron 
kinetic energy. On the other hand, at 2. 5 MeV the Born approximation gave 
relatively good agreement with experimental results. 

Screening effects due to the atomic electrons were found to be 
insignificant in the energy range of interest, except for photons with a low 
percentage of the incident electron kinetic energy. Therefore these effects 
were not included in these calculations. 

A correction factor to the Born approximation was derived by 
comparing the two theoretical approaches. This factor is a function of the 
photon share of the incident electron energy (k/T ) , the velocity of the 

incident electron (P q /E ) , and the atomic number of the scattering media. 

The corrected Born approximation was then calculated for incident 
electron energies from 0.2 to 2.5 MeV, and for a range of atomic numbers. 

It was found that the corrected cross sections gave relatively accurate 
results at the lower portion of the energy region under investigation; while, 
at energies above this range, e.g. , at 2.5 MeV, the uncorrected cross 
sections were more accurate. 

The correction factor failed to give accurate results after the Born 
approximation passed through the "transition region. " This was most likely 
a result of the limitation of the data available from the partial wave method. 
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At the present time, this method has an upper limit of around 1. 0 MeV, 
However, should improvements be made in the computer program, it may be 
possible to obtain a correction factor that is more accurate in the higher 
energy range than the present case seems to indicate. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, April 27, 1970 
124-09-11-14 
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APPENDIX A 


ELWERT COULOMB CORRECTION FACTOR FOR LOW ENERGIES 

To correct for the Coulomb effects not included in the born 
approximation, Elwert [ 9] developed a correction factor valid up to an 
energy on the order of 0. 1 MeV. This factor was derived by applying the 
techniques of mathematical analysis to a comparison between the non- 
relativistic Born approximation and the exact Sommerfeld results. It is 
therefore a semi-empirical factor limited to electrons with velocities in the 
nonrelativistic range. The Elwert factor may be written as 

/3 q { 1-exp [-(2 ttZ/137 j 8 o )]} 
f E = j3 {i-exp[ - ( 27 tZ/ 137/3)] ) 


The corrected cross section is then obtained by 


l 

\ dk /corrected 



Born 


Figure A-i(a) shows the magnitude of this factor for a range of 
atomic numbers, photon energies, and incident electron kinetic energy. 

Figure A -2 illustrates the effect of the Elwert factor on the Bethe- 
Heitler differential cross sections and makes a comparison with the experi- 
mental results obtained from Reference 7. 

Note that f — i as k — 0. This is true since f £ approaches 1 
asymptotically as /3 — < as k — ► 0. 

The maximum energy that can be radiated in a bremsstrahlung 

collision is limited to the incident kinetic energy, (hv) = T . This is 

max o 

known as the "high frequency (short wavelength) limit. TT As this limit is 
approached, the Bethe-Heitler theory breaks down since, contrary to 
experimental data, it predicts a vanishing cross section [5] . The Elwert 
factor does not correct for this. To illustrate this, one first notes that the 
Elwert factor [5] is valid only if 
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k = 0.7T 
(Z = 79) 


\ k ■ °- 5 T» 


k = OJT 


INCIDENT ELECTRON KINETIC ENERGY (MeV) 

Figure A-i(a) . Elwert Coulomb correction factor — effect of atomic 
number and photon energy on magnitude of correction. 


HIGH FREQUENCY LIMIT 


k/T o (%) 


Figure A-i(b) . Elwert factor for 0. 3 MeV incident electrons 
as a function of photon fractional energy. 





(Z/137) (/3 _1 -/3 o _1 ) « 1 . 


Since p = (E 2 - i) 1 / 2 and p = (E 2 - 1) , one may write 


-1 E 


E 


2-1 \ 1/2 


and jS 


-l 


E 


(EM) 


(E 2 -i) V 2 
o 


then our restriction becomes 


Z/137 


9 -l/ 2 

E(E 2 -1) - E q (E 


2 _ 1 ) 1/2 
o 


« 1 


Now E^ = k + E by conservation of energy; and, since E q = T q + 1, 
then E = T q + i - k. As the high frequency limit is approached, k — T q , 
then E — 1 and one has in the 


lim 
k — T 

c 

E — ► 1 


Z/137 


- 1/2 , -i /2 

E(i - 1) - E (E 2 - 1) 

0 0 


« 1 


Thus, the second term approaches infinity and the Elwert factor is 
not valid at the high frequency limit (Fig. A -1(b)) . 

It can be concluded from the data presented and the above argument 
that the Elwert factor is valid for lower energies and increases as the high 
frequency limit is approached or the atomic number of the scattering 
medium increases. Furthermore, the factor cannot be used at the high 
frequency limit, since it increases without bound at that point. 
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APPENDIX B 

NORMALIZATION OF THE RADIAL WAVE FUNCTIONS 


The radial functions for the continuum states are normalized on the 
energy scale. This requires that [17] 

/ d^ ^ E , + i p E = 6(E - E') . 


First, consider the asymptotic form of the radial functions. As 
r — ► oo , 


rf — -A(E - l) 1 / 2 sin (pr - 6) 

(B-i) 

rg — ► A ( E + i ) V 2 cos ( pr + 6 ) 

The asymptotic form is also given by letting r — °o in equation (C-ll) , 
Appendix C; i. e. , 


rf — i(E-l) 1 / 2 N T (2y+l) (2pr) 7 


(r + iy) exp(ipr+j77) 
T (y+l+iy) 


rg — (E+l) l/2 N T (2y+l) (2pr) 7 


x (2ipr) iy " y - C.C. 


(y+iy) exp( ipr +ir? ) 
T (y+l+iy) 


(B-2) 


x (2ipr) iy-7 + C. C. 


Now [17], 


y+iy _ exp[-i ARG T (y+iy)] 

r (y+l+iy) I r ( y +iy ) I 
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I 



and since 


( 2ipr) iy = e 


TT 

-~ y 


iy In ( 2pr) 


and 


(i)' 7 =e 


. 7 r 

1 ~o~ r 


the first term in square brackets of equation (B-2a) can then be written as 

TT 

ipr - 
e 


(2pr) 


-y 


I r (y +iy) 


exp i( -ARG T (y +iy) +y In 2pr - y + 17) 


ipr-— y 


Define 6 = y In ( 2pr) - ARG T (y +iy) + V - — y ; then, 

rf — i (E-l) 1 / 2 N r(2y+i) (2pr) y (2pr) jexp^ i^pr - y + 6^ J 

- exp £-i (pr - ~ y + 


or 


7 r 


rf 


2(E-l) 1 / 2 NT ( 2y +1) e 
I r ( y +iy ) I 


sin (pr + 6 ) 


(B-2b) 


A comparison of equations (B-i) and (B-2b) indicates that these equations 
will be equal as r — °° , if 


_ 7ry 

2 N T(2y +1) e 2 
|T(y+iy) ! 


(B-3) 


Now use the rule for normalization of the eigenfunctions of the 
continuous energy spectrum [21] 
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fp, (r) dp’ = 1 


( B-4) 


oo p+A p 

J r 2 dr f p ( r) / 

0 p-A p 


where Ap is a small momentum interval. The normalization factor will 
first be found on the momentum scale and then transformed to the energy 
scale. 


The asymptotic form of f is given by equation ( B-l) as 
f = A/r [cos(pr + 6)] 


Consider the second integral in equation (B-4) 
p+Ap 


J — cos(p'r + 6) = — cos(pr + <5) 
p-Ap 1 1 


sin Apr 
r 


where values on the order of 1/pr and A p/p have been neglected. 
Substitution of this result into equation (B-4) gives, after replacing the 
oscillating function cos 2 (pr + 6) by its average value § , 


A 2 f dr 

0 


1 . 


This integral can be found in any standard table of integrals; then, 

A ! (f) ’ 1 ° r A - J4- . 

This is the normalization factor on the energy scale. 

If f are eigenfunctions normalized on the energy scale, and f 

JCj i . 

are eigenfunctions normalized on the momentum scale, they are related by [21] 
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Thus, on the energy scale, one 


Since E = |p 2 , then 



has 



This value of A must be set equal to the value of A 


found in equation (B-3) to determine N, the radial function normalization 
factor. Hence, 


jl 

2 N r (2-y +1) 2 /_2_\ 

I T ( 7 +iy ) I 6 \7rp/ 


7ry 

N = 1 r ( -y +iy ) 1 e 2 

42¥p r (2y+l) 

, Z e 2 E 

where y = 


(B-5) 
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APPENDIX C 


THE DIRAC EQUATION IN POLAR COORDINATES 


The bremsstrahlung problem is described by a relativistic electron 
in a central field. The equation of motion for this particle is given by the 
Dirac equation using the Hamiltonian: 


H= a • p +JJ_+V(r) , 


and 

Wp = E ; p (C-l) 

[( g * ~9+&) + V(r)]| = E| , 
where a and £ are the Dirac matrices defined below and 



Since a central field has been assumed, the Dirac equation will be 
separated into polar coordinates using the method presented by Rose [19] 
(In this appendix, extensive use will be made of two books by M. E. Rose 
[17,19].) 

First one writes the gradient using the identity 
V = r(r*V)-rx(rxV) , 


where r is the position vector for the Dirac particle relative to the field 
source. Since L = -irx V , where L is the orbital angular momentum 
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HMI 


operator, one has 


■=* — 9 . r -* 

V = r 1 — x L 

9r r 


( C-2) 


For any two vectors A and B whose components commute with the 
components of g , the 2x2 Pauli spin matrix, one has 
(<T * A) (a - B) = A* • B - + i cT • (AxB) . 


Now, a = 



and £ = 


4x 4 


' I 0 


0 -I 


4x 4 


where each element is a 2 x 2 matrix. One now makes the following 
definitions: 


I = the 4x4 unit matrix 



= 0 

= the 4x4 matrix ( J , where each a is a 2 x 2 matrix. 

0 <T 

\ “ 

For example, 


u- 



, etc. 


P 3 = 
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Note that 


? = E i ?3 = £3 E and £3 P3 = I 


-1 

which implies that p 3 = _p 3 . 


Consider g • A g • B , where A and 
components commute with the components of g 


B are any vectors whose 
One has 


aT- A o’- B =£3 E ' ^*£3 E ' 

= £3 £3 E • r E • b • 

Now, since p 3 = p 3 1 , then 


O' 


A o; • B = 





where the subscript 4 indicates a 4x 4 matrix. 
However,- 


ct'Act' B = A • B + i g • (AxB) ; 


then, 

a • A a • B 



+ i 


g • (A x B) 0 

0 g • (A x B)/ 4 


Let A = r and B = L ; then, by using a vector identity 


r • L = r • (-irx V) = [rx( -ir) ] • V = 0 


one has A • B = 0 . 
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Thus, a • r a • L = i 


g • ( r x L) 0 


0 g • (r X L)/ 4 


0 g • (rx L)\ _ 

= iPs | _ _ _ ) =i ^3 J 

g • (r x L) 0 / 4 


and a • r p 3 a • L = ia ■ rx L , since a p 3 = p 3 a . 

Therefore, a - r ^*L=iQ!*rxL. Now, ^ • r = a I 

where a is a 4 x 4 matrix and ia • rxL = ff ) ,• L ; then 
= r = = r y 


i a • r x L = a 
= = r 


<j • L 0 


0 a * L, 


^ ^ g ^ ^ 

Also, since a • V = “ “ £1 ' r * L (from equation C-2) 


then a; • V = o; — - a 
= =r 3r r 


a • L 0 


0 cr • L, 


and 


a; • V = a; ( — Z ' L 

= =r \ 8r r “ 


The Dirac equation can now be written as 


(E-V)j £ = 


>(£-T + 


Now, define the operator [19] 


Ks S(1 + Z ' L) , 


(rx L) 


= a , 
- r 


( C-3) 


( C-4) 
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or 

r-(£ • r ) +2 (f • r )*i- ( r+ I /2 ) + T i • 

Now, ^l7+ Ij/zJ 2 = J 2 , the total angular momentum operator, with 
eigenvalues j(j + 1 ) . 

Define the eigenvalues of K to be -k ; i.e., 

K 2 J = k 2 | , 

or 

+ j^ 2 + | = j(j + i) + -j- $ = (j +-|-y £ . 

Therefore, the eigenvalues of K 2 are /c 2 j + 

k = ± i, ± 2, ± 3, • • • or j = |k | - - 5 - , k =£ 0 ; where, for the upper 
spinor (equation C-5) 

f = K if K > 0 

i = -K-i if K < 0 . 

This could be written as l = j + — - S(/c) , where S(k) is the sign of 
kappa; i.e., S (k) =k/\k\ . 

To indicate the dependence of j and i on k, the eigenvalues of j 
by j and I by 1 in the upper spinor are replaced. For the lower spinor 

K K 

(equation C-5) , define 
7 =i 
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where 


i 

= K “ 1 

for 

K > 0 

“1C 



i 

= -fC 

for 

K < 0 


~K 


That is. 


i = l ~ S(»c) . 

~K K 


The wave function for the Dirac equation is written as 



where £ and $ are the two-component, single-column matrices referred 
to as the "upper" and "lower" functions, respectively. The wave functions 
are proportional to the angular functions, ft and ft , respectively, 

where the eigenvalue ju is the z component of the total angular momentum. 

The angular functions may be written as a product of the orbital 
angular momentum function and the intrinsic spin function, coupled by the 
vector addition or Clebsch-Gordon coefficients: 


ft 

} ft 


Z 



ju-m, m 


Y M-m 

i 

±K 


A 


(r) 


2 * 


m 


where the two-component Pauli spinors are defined as rj x x 
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u 

One expresses the wave function ^ as 


if- 


g(r) S2 
if( r) £2 


KH 


-KM, 


4x i 


( C-5) 


where f and g are the radial wave functions . 

Substituting ip into equation ( C -3) one has, after noting that 




where the primes denote differentation with respect to r. 


Since g r - 3 and g r (=^ ' gr (») 4><1 


the Dirac equation becomes 


'g S2 


(E-V) 


(CJ u 


if £2 


Z 


'if' £2 


-K/i\ 


<=! r 


-KJUy 


g’ S2 


4x i 


KJU 


4x 1 


Sr / if J-LS2_ 


+ i 


KM\ 


g S ’ L £2 


KM 


g ft 


KM 


-if £2 


4x 1 


-KMy 


4x 1 
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From equation ( C-4) , 


K = £ ( i + 2 • L ) = 


'(1 + 2 ‘ L ) 


-(1 + a • L) 


which implies that 
£2 


K 


~KMl 


£2 


KM 


4x 1 


K £2 


“KMl 


[-K £2 


KM, 


4x 1 


since K ip = - k . 
Then, 


(1 

+ 

Of • 

L) 

£2 

= K 

£2 



-KM 


-KM 

(1 

+ 

Of • 

L) 

£2 

“K 

£2 



KM 


KM 

One 

next considers 

a 

£2 






= r 

KM 


to find the results of carrying out this operation. Now a is a pseudoscalar 

operator so that j 2 commutes with cr . That is, for i , 

= r J z 


r(j , a ) = (i +S , 
J z =r z z 


x a + y cr ) = 0 
= x J =y 


since 


(i , x) = iy, (i , y) = -ix, (S , a )=ia , 
z z z — x — y 


and (S , cr ) = -i cr . Now, since cr has odd parity, 
z =y =x =r ^ ,r ’ 
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tr = a £2 , 

= r K/Lt -KfJL 


where a 2 = 1 , since <7 2 = 1 . The latter follows from the result 
=r 

a • Aff ■ B=A *B + iff(Ax B) , which implies that 

o;^ = a • rg • r = r • r + 1 a • (rxr) = l z + ia(0) =1 . Now take r along 
the Z axis; then 0=0. The spherical harmonic in the definition of £2 . 
becomes 


m 


A 

(e 3 ) = 


21 + 1 

K 

, 47T 


1/2 


m, 0 


Setting m = 0 will then give 

(21 + t\ V 2 ^ j K 

£2 = l ^ ) C A V 

•1 jj l g j 

Therefore, a(2f +1) 1 / 2 C . K K -2y(2{. +1) 1 / 2 C. K ,v 

K Op K 0 JLl 

four possible cases: 


There are 


j = i ± 5 , M =± 5 . 

KK 

One finds, after substitution of these values and solving for a, that a = -i 
for all four possibilities. Therefore, 


cr £2 = - 0 

= r KjU -k/u 


and, by the same type of argument, one can show that 

u 12 = - £2 

= r -icy icy 
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After one operates on the angular functions with the operators 
(1 + a • L) and o - , one's set of equations becomes 

[ (E-V-i) g + f' + f/r - k f/r] Q =0 
i [ ( E -V+i) f - g* - g/r - /c g/r] fi _ =0 . 

~K fJj 


This implies that the expression in square brackets must equal zero; i.e. , 


(E-V-l) + f' + f/r - k f/r = 0 


(E-V+l) - g’ - g/r - k g/r = 0 


(C-6) 


For simplicity, one writes Uj = rg and u 2 = rf; then equation (C-6) 
becomes 



- u A + (E+i-V) u 2 


du 2 

dr 


u 2 - (E-i-V) Uj 


(C-7) 


These are the radial equations of a Dirac particle in a central field. The 
boundary conditions required are that the resulting solutions must remain 
finite as r — 0 and as r — . 

The asymptotic forms of equation (C-7) are found by assuming that 
the potential tends to zero sufficiently fast for large r (i.e. , V 0 as 
r — °o ) . 

For large r, equation ( C-7) becomes 



du 9 

dr 


+ (E - 1) u t = 0 . 
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The general solutions for these equations are 


=N(e" Xr + q(E) e Xr ) 

and 

„ T /1-eV/ 2 -Xr /1-EX 1 / 2 

U2 = N [’\JTe 7 6 +q(E) (W ) 6 

where X = (i-E 2 ) 1 / 2 , 


(C-8) 


N 

q(E) 

Since 
then becomes 

i Ul = 

J 

u 2 = i 


is the normalization factor, and 

is found from the condition that the solutions remain finite 
as r -*■ 0. 

E > 1, let p = (E 2 - i) 1 / 2 so that X = ip ; equation (C-8) 
i^ N ( E + 1 ) 1//z sin ^pr “ 4 ^ + 

N(E - i) 1 / 2 sin ^pr - + 6 ^ , 


where p is the momentum of the particle and 5 is the phase shift that will 
depend upon the nature of the external field of the particle. 

Following the method of Rose [17], one lets V = -ze/r = Z/r in 
equation (C-7) ; i.e., 
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After the substitutions 

ut= (E +1) 1 / 2 (<t> i + <p 2 ) , 

u 2 = i (E - l) 1 / 2 (4>i - 0 2 ) , 

P - (E 2 - l) 1 / 2 , 

and 


x = 2ipr (where (j) i and (f> 2 are functions of x) are made, 
one then has 


= /i_ +i Z E \ /jc 
dx \ 2 p x \x 



4>2 


and 



After taking the complex conjugate of these equations, i.e. , 



and 


d0 ? * = _ /j£_ _ 
dx \ x 




+ 


iZE 

px 


4>2 




one can see by inspection that the two sets of equations are equal if 
t p * = <p 2 . After eliminating <p 2 , one has 



where y = [/c 2 — Z 2 ] 1 / 2 . 
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Let M = x 2 <pi to get equation ( C-9) into the standard form of 
Whittaker’s equation 



Equation (C-9) becomes 


d z M 

“d?" 


i_ + (± + 1 ™)± 

4 V 2 p / x 


y - 


X 


= 0 


This equation has the regular solution at r = 0 : 

,/y+i -x/2 


M(x) — X ' 


F(y + i + iy, 2y + 1, x) , 


where y = ZE/p and F(y + 1 + iy, 2y + 1, x) are the confluent 
hypergeometric functions [22] 


> = tr + - 


oo a m 

E m _x 

a c m - 

m=0 m 


Now, M = x 2 <p t ; then 

(p 1 =^x y e X// 2 F(y + 1 + iy, 2y + i, x) 

= N (y + iy) e 1 ^ (2pr) y e l ^ V F(y + 1 + iy, 2y + i, 2ipr) 
= N(y + iy) e 1T? (2p) Y 0(r) , 


where N is the normalization factor. Next, the phase of rj is chosen such 
that <p 2 = ; this requires that [17] 
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- 2iT ? = _ Y + i y £ 1 L ^ i 

7 - iy k - iy/E 0* dr \ pr ) cp* 

This equation can be evaluated by using Kummer's formula 

e X//2 F(y + i + iy, 2y + 1, x) = e X ^ 2 F<y - iy, 2y + 1, -x) 

and 

xF(a + 1, c + 1, x) = C [F(a + 1, c, x) - F(a,c,x)] 


One finds that 


2 i r] _ k - iy/E 

y + iy 


(C-iO) 


The radial wave functions may then be written 


rf = i(E - i ) 1 / 2 N(2pr)^ {(y + iy) exp( -ipr + i tj) 


rg = (E + i ) 1 / 2 N(2pr)^ {(y + iy) exp(-ipr i tj) 


x F(y+l+iy, 2y+l, 2 ipr) - complex conjugate} 

(C-li) 


x F(y+l+iy, 2y+l, 2 ipr) + complex conjugate) , 


where N is the normalization factor and rj is the phase determined from 
equation (C-10). 

The solutions are normalized in the energy scale. If ip and ip 

-b ih ' 

are solutions corresponding to energies E and E', then one requires that 
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/ d 3 x^ t 4 > e = 6(E -E») . 


From Appendix B, equation (B-5) , the normalization constant is 

7ry 

N = lr <^y - )l e 2 

2(irp)‘/ 2 r(2y+l) 


Therefore, the radial wave functions for a particle in a Coulomb field 
normalized on the energy scale, and with energy corresponding to the 
continuous energy spectrum, are 


(C-12) 


II 

y 2 


rf = 


rg = 


i(E - l)*/ 2 (2pr) 7 e 
2(7Tp) 1 / 2 r ( 2y+i ) 


vry 

( E + l) 1 / 2 (2pr) y e 2 
2(7r p ) 1 / 2 r (2y+l) 


I r (y + iy)| {exp( -ipr + i tj) 


x (y + iy) F(y + 1 + iy, 2y + 1, 2ipr) - C.C. } 


I r (y + iy) I {exp( -ipr + i i?) 


x (y + iy) F(y + 1 + iy, 2y + 1, 2ipr) + C.C.} 
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